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INTRODUCTION 


Every beginning course In classical mechanics has as an 
exercise the derivation of the Rutherford scattering cross 
section for two classical Coulomb point charges. Many people 
are surprised to learn that this is still apparently an un- 
solved problem if a uniform dc magnetic field is assumed to 
be present. Since the interaction of classical point charges 
is fundamental for all discrete-particle processes that go on 
in most plasmas, this simple example should make clear that 
there are many things we do not know about such processes. 

The situation becomes worse if we add in the complexities 
associated with the long range of the Coulomb force, and the 
statistical mechanical side of the problem. 

There is a seriously underpopulated border region of plasma 
physics that lies between controlled thermonuclear research 
and the older, less volatile society of classical statistical 
mechanics. This territory was initially explored in classic 
papers of Landau [1,2], Vlasov, [3], and Bogolyubov [4], Ten 
or twelve years ago a few permanent settlements appeared to 
have been established there. Life on this rugged frontier 
proved to be harsh, however, and many of these settlements have 
now been abandoned. One hope of these lectures will be to re- 
cruit some vigorous pioneers who will be willing to try again 
at the neglected task of developing this important area. 

Our subject matter may be roughly defined as that class 
of plasma processes for which the so-called "Vlasov approx- 
imation" is inadequate [5], Such phenomena include: equations 

of state and other equilibrium thermodynamic relations, thermal 
relaxation phenomena, transport properties such as diffusion 
and electrical conduction, and microscopic statistical fluc- 
tuations in such quantities as the electric field and the 
charge density. All these may be loosely called "discrete 
particle" processes. 

Convincing calculations of these quantities have been 
given in some cases for the case of a plasma in the absence 
of an external magnetic field. Calculations in the presence 
of a strong magnetic field are usually absent or less than 
convincing, though progress has been made within the last year 
or two. Reliable laboratory measurements of these quantities, 
that would pass muster in other branches of physics, are 
usually lacking, both with and without the external magnetic 
field. 
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Results In this area come slowly, and it is rather easy 
to make embarrassing conceptual mistakes, even though it is 
only classical physics we are dealing with. There is also 
great economic pressure for plasma physics to show results in 
a hurry, more rapidly than it is realistic to expect these 
rather deep fundamental problems to sort themselves out. To 
yield to this economic pressure is to gamble upon achieving 
practical success in plasma situations without needing to 
understand these phenomena. Such a gamble may be justified 
for some people, but I think it would-be unwise for alt of us 
to take it. 

The question of why practitioners of orthodox statistical 
mechanics have largely stayed away from these problems is more 
difficult to answer. Their training and extensive experience 
with similar problems in neutral particle systems would ap- 
pear to be an excellent preparation for plasmas. Perhaps it 
has to do with their perceptions of the sometimes desparate. 
professional style of the plasma physics community. But for 
whatever reason, the classical statistical mechanics experts 
have largely remained aloof from considerations of systems witt 
long range forces, one suspects much to the detriment of both 
subjects. 

The subject of discrete-particle processes is much better 
understood in the absence of an externally-imposed dc magnetic 
field than in the presence of such a field. Many, if not most, 
of the important laboratory plasmas involve external magnetic 
fields in a fundamental way. The emphasis in this article 
will be almost exclusively on the situation in which a strong 
external magnetic field is present. A recent volume by the 
author [6] surveys the theory of the unmagnetized , or field- 
free plasma, and many other good survey volumes exist. 

The following material is divided into four sections. 

In Section I, a number of results from the equilibrium sta- 
tistical mechanics of two-dimensional plasmas are derived. The 
corresponding results for three dimensions are already well 
known. These results are independent of the presence of an 
external dc magnetic field, and are of interest in their own 
right as well as being essential to what follows. Section XI 
is concerned with the non- equilibrium statistical mechanics of 
the electrostatic guiding-center plasma, a two-dimensional 
plasma model of great power and simplicity, recently Introduced 
by Taylor and McNamara [7]* Section III concerns the generali- 
zation of this model to three dimensions. Section IV returns 
to two dimensions and relaxes the guiding- center model to in- 
clude finite Larmor radius effects. 
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I. EQUILIBRIUM ENSEMBLE AVERAGES ; 

STATIONARY CORRELATION FUNCTIONS 

(A) Phase Space and the Gibbs Canonical Ensemble 

As is the case with any system with a very large number 
of degrees of freedom, our description is necessarily proba- 
bilistic. The dynamical state of the system can be specified 
by identifying it with a point of the phase space of the sys- 
tem. The coordinates of any such point (X, say) are the posi- 
tions and velocities of every particle. The dynamical state 
of the system is then equivalent to a random variable defined 
over the phase space. Equilibrium statistical mechanics [1] 
demonstrates that the normalized probability distribution of 
this random variable appropriate to a system in thermal equi- 
librium with a temperature 0 = KT (in energy units) is 

I). = exp(-s/9) 

e< l* J dX exp(-£/e) 


(i) 
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This is, of course, the canonical distribution of Gibbs; £ is 
the total energy, uniquely defined for each point X of the 
phase space, and dX is the volume element. 

Thermodynamic quantities and other macroscopic observables 
are interpretable as "ensemble averages", or expectation values 
computed with respect to D eq . Thus, if A(X) is any dynamical 
variable, its expectation value will be written as 


<A> = J dX A(X) D eq< 


( 2 ) 


where the integration runs over the entire phase space acces- 
sible to the system. 

Thermodynamic functions are derivable from the partition 
function , which, apart from quantum mechanical factors which 
will not be of interest to us here, is 

Z = J dX exp(-e/0) (3) 

In particular, the equation of state for the pressure [2] is 


P = 


av 


jU/Z 


0_ 5Z 

z sv 


(4) 


where V is the spatial volume in which the system is confined. 

If one wants to calculate extensive quantities (such as the 
Helmholtz free energy or the entropy) accurately, one has to 
put the additional factors into Eq. (3) correctly, but since 
they are volume-independent, they cancel out of Eq. (4). 

For definiteness, we will usually consider N positive 
charges +e and N negative ones -e. For e we have the expression 


£ = T + U 


(5) 


T = 




( 6 ) 




( 7 ) 
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T is the total kinetic energy and U is the total potential 
(interaction) energy. The^I^ in (6) runs over all particles 
of both signs, and m± and are the mass and velocity of the 
ith particle. In (7), is the two-body Coulomb potential 
between charged particles i and j . The sum runs over all pairs 
We make the convenient and important restriction to electro- 
static interactions , ignoring all retardation effects and mag- 
netic interactions as being down by a factor of 0(v 2 /c 2 ) from 
the ones we shall include. (This approximation is often called 
somewhat incorrectly, the "low beta" approximation by plasma 
physicists. "Beta" is a rather important plasma quantity which 
is essentially the ratio of thermal energy density to magnetic 
energy density.) It is also important to note that the things 
that have been said so far are unchanged by the addition of an 
external, constant, and spatially uniform magnetic field of 
arbitrary strength, since & is unchanged by the addition of 
such a field. Thus, the calculation of such ensemble averages 
as Eq. (2) are independent of l£, the external field, so long 
as A is a time-independent function of X. 

(B) Reduced Probability Distributions [2-6] 

It is convenient to project the probability distribution 
(1) onto a subspace which is just the phase of s particles, 
where s is some small integer (1,2,3,...). Let X * E (x,,v.) 
be the position and velocity of the jth charge. Then J J 





r 


dX 




D 

eq. 


( 8 ) 


is the reduced probability distribution in the phase space of 
s particles. The volume factor V s is for later convenience, 
and the integration is over all the coordinates of all the 
remaining charges. To discuss two charged species, it is con- 
venient to partition s into s-^ "ions" (positive charges) and 
s e "electrons" (s e + s.j_ = s) and define 
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»eq. 

1 





= V 


s.+s A 
i e 


J 


dX 


a +1 




dX 


8 e +1 




eg. 


where X 




(9) 


,+i ->-i, 

‘vV 


are the phase space coordinates of the 


jth ion; a similar expression applies for electrons. The 
superscripts will indicate the species of the charged par- 
ticles. Note that all the f are symmetric under the 

3 i ,S e 

interchange of like particle coordinates. 


It is immediately obvious that the velocity-space parts 
of the f are trivial and may be separated from the con- 

S . a S 

i e 

figuration-space parts : 



6 


e 



expC-ity?^ /2e)"'l 

(aire/m^ 2 J 



exp(-m e v: /2e)"'l 

W ( 


(2TT0/m e ) 


J 


s. ,s 

e 


( 10 ) 


where d is the dimensionality of the system, and 
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n a 

V 6 e 


1 e / d5 s.+i *•* +1 d3 S ex p(- u / 0 ) 

^ 1 6 

J dx^ ... d5c^ dSc^ ... dj^ exp(-U/6) 


( 11 ) 

is the configuration-space probability distribution of the sub- 
system. 

It will be generally understood that we are dealing through- 
out with very large systems: N » 1, and N/V = n Q , the number 

density. The "thermodynamic limit" of N «> and V 00 may or 
may not be well defined for such expressions as Eqs. (3) and 
(4) . It is not an entirely academic question as to whether such 
limits exist. Some rigorous work has been done in one dimen- 
sion [7,8] and in three dimensions [9], and more is badly needed 
(in two dimensions, especially). We shall generally assume that 
the limit does exist; but that V and N are just very large but 
finite quantities is perhaps a slighly more acceptable way for 
us to look at the problem. The large-system limit is desirable 
in that we may hope that removing the system's boundaries to 
infinity will leave us with a "bulk" limit in which the surface 
contributions to extensive quantities will be negligible com- 
pared to the volume contributions, and that nothing that hap- 
pens near the boundaries will affect the intensive quantities. 
(This may appear to be obvious, but it is not.) Making this 
assumption, it is easy to show that n 1 Q - n 0 l - const. - 1, 
and that the higher n depend only on the separations of 

V 3 e 

the particle coordinates and not on their absolute spatial lo- 
cations . 

(C) The Equilibrium BBGKY Hierarchy 

Most of the macroscopic functions of interest can be ex- 
pressed as integrals of n 2 Q ,n Q 2 , and n x x , as it turns out. 

Any approximate method of calculating these quantities is 
therefore of some interest. A convenient starting place is 
the well-knwon equilibrium BBGKY [2-5] hierarchy , which can 
be readily obtained by taking gradients of Eq . (11). It is 
left as an exercise to show that, assuming periodic boundary 
conditions over a very large volume V, for any j between 1 
and s . , 

l 
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S +3 
e 


The notation used is the following : <j> is the Coulomb poten- 
tial interaction of the twt> ions, and (|) ee is that of two elec- 
trons. (pKx) = cf> ee (x) . = <() e ^- — — $ ee is the Coulomb po- 
tential between an electron and an ion. ^(s^s^) is the total 

Coulomb potential energy of the first ions and s e electrons 

If s e and si are not large, we can replace (N - si)/V and 
(N - s e )/V by the average number density n Q , on the right-hand 
side of (12) , A similar relation comes from taking gradients 
of Eq. (11) with respect to the electron coordinates: 



ai(s i ,s e ) 
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r 



,-*e 

dx 


s +1 


Sep 


ee 




-*e 

X s +1 
e 





a +1 
e 


(13) 


where we have assumed N » or s e - The normalization con- 
ditions to be obeyed by n are, from (11) , 

S - 9 s 



dx 


= V 


s . +s 
x e 


(14) 


and n o , o n i ,0 


= n 


o,i 


= 1 . 


(D) Expansion in the ft Plasma Parameter" e 

For the entirely general case, there is no particular ad- 
vantage in writing the thermal equilibrium relations in the 
form of Eqs . (12) and (13). Rather, their advantage lies in 
the ease with which perturbation theory can be done on them. 

It turns out that for a wide range of parameters, the poten- 
tial energy of two charges in a plasma which lie at a distance 
equal to the nearest-neighbor separation is much less than the 
average kinetic energy per particle, 6. This means that the 
Interactions are mostly weak in the sense that 4> i: V0 

and cf> ie /6 can all be considered to be « 1 over most of the 
phase space. We can thus do a weak-coupling expansion in 
<< J )> avg.^» (where <4> > aV g. is some average interaction poten- 
tial) . The expansion is complicated slightly by the long range 
of the potential, which raises some of the integral terms in 

Eqs. (12), (13) by one order in <<{>> /0. 

avg ■ . 

The details of this expansion for the three-dimensional 
case have been adequately described several times (see e.g. , 
Guernsey [10] or Montgomery [11]) , and there is no purpose in 
doing so again here. We shall outline the development, how- 
ever, for the two-dimensional [12] case, since so much of the 
later theory to be developed involves two-dimensional models. 
Equations (12) and (13) apply to one, two, or three dimensions, 
with appropriate modifications of the Coulomb potential 
<f>ii ss (j>ee = ^ei = -(j>ie, in tw o dimensions, the 'charges are 
imagined as very long parallel "rods" of length £ and charge 
e, and the potential between two rods of positive charge +e 
each is, for example, 
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ii,-*i ->iv 

q> (x^ - *j) - 


2e 


„ i *4i -*i 

** i *« - x j 


(15) 


with similar expressions for the other interactions. 

The zeroth order (i.e., non-interacting) values of the 
n are seen to be those solutions of Eqs. (12) and (13) 

s i ,s e 

with all the potential energy parts neglected and which also 
satisfy (14) . These are, obviously, 



s 

e 


= i, 


all s. and s 
x e 


(16) 


The superscript zero means of zeroth order in powers of the 
interaction potential, and we imagine (16) as the first term 
of the series. 


V,S =n i°]s 
l * e a/ e 


+ n 


(i) 


(17) 


We expect the higher-order terms of this series, the Correla- 
tion functions 0 , all to vanish as the particle separations 
become very large, and n n(0) then. 

s i> s e s i’ s e 

The lowest order there is, according to (16), no correla- 
tion between the positions of the particles: what one would 

expect from non-interacting particles. We proceed to (12) 
and (13) to obtain the next terms in the series. However, the 
long range of the Coulomb potential means that one has to keep 
the next order terms in the integrals in Eq. (12) , above those 
a purely formal estimate would indicate. 

A self-consistent expansion procedure is readily found 
in which the correlations between larger numbers of particles 
is progressively weaker in powers of <<|)>/0. The next order can 
be found in terms of pair correlations: 



MAGNETIZED PLASMA MODELS 


441 


+ ) ee<-+e ne> 

L . • P (x, - x ) 

k<4=l K 1 
b 4 s 

♦I 

k=l 4=1 


(18) 


where 


ii 


P 


ee ei 
P = ~P * 


With the assumptions, it can 
part of all of Eqs. (12) and (13) 


p ee (x) = e 2 £ftx) , 0 ee (x) = e 2 ^(x) , 


be shown that the first-order 
is satisfied if one sets 

and 




a x, 


12 


1 

6 


a +(* 12 ) 

*♦ 

a x 12 



0 



a K x 13 ) 

a x i? 


^(x 25 ) 


(19) 




a H x 13 ) 

5 x l3 


<52( x 23 ) 


• * y I 

In Eq. (19), x 12 = |xi - x 2 | , and so on. It is a differentio- 
integral equation of the convolution type, and can be solved 
by the Fourier transformation. Setting 
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$ (x) s I die (P(x) e lk ' x 


ty(x) = dS *(g) e ilC X 


<9(t) = 


2n e 

1 + | — ( 2tt ) 2 |( k ) 


P ee (3?) 


^ e (g)/9 . 

l+(2 n o /0) cp ee (lt)(2 n) 2 


Poisson’s equation, Fourier-transformed in two dimensions, give 

4> ee (k) = e 2 /Tr£k 2 , so that we may finally write (22) as the 
classical Debye-Hiickel expression : 


p ee (£) = 


8 tt 2 ii. 


- , 2. 2 
1 + k \ 


Here, we have introduced the very important Debye Zength into 
the problem, namely 


M = 

2 ” 


8 Tin e 
o 


If we Fourier transform Eq, (23) back to x-space, we get 


p (x) = - 


K 0 (*/V 


(25) 
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where K 0 is a Bessel function of imaginary argument, and has 
the asymptotic form 

P ee <*> — -> - *£■ 7 -sr <”® <-*/V <26) 

x » v 


We see that X D has the important physical significance of being 
the distance over which pair correlations persist, just as in 
three dimensions, though the functional form is somewhat dif- 
ii ee ei 

ferent. Since p = p = -p , we now have all the first- 
order solution* 


The calculation of the correlation functions to higher 
order becomes quite messy, but as long as p ee (x) is << 1, Eqs. 
(17), (18), and (25) are a satisfactory approximation* It Is 
readily seen that the condition for validity is thus 2e 2 /£0«l, 
or what is equivalent, 

? - 1 

e ■ (4nn « 1 (27) 


This condition is essentially that the number of charged. rods 
per Debye square is very large. The well-known analogue [13] 
of this condition in three dimensions is that the number of 
point charges per Debye cube is very large. 


(E) Further Comments on the Validity of the Expansion 

The failure of the weak- interaction approximation over 

small regions of the phase space (at small particle separations) 

is indicated by the fact that the short-range form of Eq. (25), 

p (x « A ) * -( 2e 2 /£0) ln(x/A_) . This ceases to be « 1 at a 
ee v u 

distance x given by 


2 e 




mm. 


= * 3 > 




(28) 


Compared with the nearest-neighbor separation n 0 , 


this gives 
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^ exp (-1/0 (29) 

n * 
o 

which is smaller than any power of e , but is not zero. Be- 
cause the separation of two charges is so seldom less than 
x min. » it is reasonable to hope that the important features 
of the theory will be unchanged by ignoring the short range 
of <j> ee (x) . To render certain integrals convergent, it is 
often convenient to work not with (f> eG (x) , but with another 
function which differs from it only for x < » namely 


Jo 



e 2 e ikx cos 0 

TT Jjk 2 


(30) 


where k max = l/Xj^ . The length x min has the physical 
interpretation in the unmagnetized case of being the distance 
of closest approach of two colliding thermal particles with 
zero impact parameter. To remove this somewhat unsatisfac- 
tory state of affairs at short range requires quantum mech- 
anical modifications in three dimensions, and has not been 
carried out in two dimensions. It has been concluded by every- 
one who has seriously looked at the three-dimensional calcula- 
tions that the dominant effects are correctly obtained from 
the modification of <f> ee of Eq ♦ (3) at short distances. We 
shall have to hope that the same thing is true in two dimen- 
sions. The effect of this modification of the interaction on 
the pair correlations is just to^cut off the Fourier integrals 
of such expressions as (23) at |k| - l/ x m i n 

(F) Ensemble Averages of Macroscopic Functions [14] 

To summarize our results to this point, we have calculated 
the four pair correlations in this problem in terms of . a single 
function, which in Fourier transform language reads 

iO/r^\ 

P W = P (k) = -p (k) = -p (k) 


(8 n 2 n o ) -1 

— ■ 

l + k S 


( 31 ) 
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This is not an exact result, and is valid to the extent 
that the inequality (27) holds. We have immediately from it 
expressions for all the various , from (17) and (18), 

s i jS e 

to the same accuracy. We can calculate, in terms of it, ad- 
tional quantities of physical Interest. For example, the 
electric field 


- r + 
hu I x “ 


N , ^ 

E (¥ ^ I 3 - 3 ? ij > 

k=l I 


can be readily shown to have zero ensemble average 


<3(x)> - j D^ 2 dX = 0 


as can the net electrical charge density 

p(x) = Y ]T (3c - j£) - 6(x - JJ)] _ (34) 

0=1 

For, using our previous expressions for p ee , p 1 , p ie , we can 


also show 


<p(x)> = 0 


However, the vanishing of the ensemble averages of E and p 
should not^be construed to mean they are zero, because, for 
example, <E 2 > ^ 0. Let us calculate the autocorrelation of 
the charge density, 


<p(x) p(x + r)> = J* D eq. P<*) p(* + r) dX 
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For this purpose it is convenient to insert (35) , (9) , and (10) 
into (36) , getting : 

i 

<p(x) p(x + r)> 



Using ni^o = n 0jl = 1, and passing to the limit of N,V very 
large, we use (18) and our previously derived expressions for 
the pair correlations to write: 

2 e 2 2 n e 2 

<p(3?) p(5f + ?)) = 4 J%- p ee (?) + — 6(r) 08) 

l l 

which depends only upon the separation of the points r, as we 
know it should. 
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It is often useful to decompose expressions like Eq. (38) 
into a Fourier sum or integral, and write them in. terras of 
their spectral densities . Thus, 


<p(x) p(x + r)> = J" die S p (k) e 1 ^ r 


(39) 


where Fourier transformation ^f Eq, (38) and use of (23) shows 
that the spectral density S^(k) is 


v*> 


n o e 

_ 2 „2 

2tt l 


r iA ] 

v iik2 4J 


(40) 


Of considerable interest also will be the spectral j|en^ity of 
the auto-correlation tensor of the electric field, S E (k) , de- 
fined by 


< 3 ( 2 ) tit 


?)> 


-/ 


dS 




lit 


(41) 


Two applications of Poisson’s equation shows that S and S 
are related by 


S E (k) - 


1 6 TT^ lilt 


s 0 (S) 


so that 


v*> 


J&7T 


Kg 

2 

k 


2 2 

1 + k 


(42) 


(43) 


It will be instructive for the reader to compute other auto- 
correlations as exercices, such as the current-current auto- 
correlation, <J(x) J(£ + r)> 
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(G) Field Energies, Thermal Energies , and a Warning 


It is interesting to calculate the rms electric field 
strength 

r r w ■ 


<2 2 > = Trace dS*? B (2) = ~ 


k dk 


, . 2.2 

1 + k 


(44) 


Equation (44) contains a logarithmic divergence at large k, 
which is to be associated with small particle separations; it 
can be rendered finite by the cutoff (30). However, this is - 
misleading, for (44) contains in it a term proportional to 
the infinite electrostatic self-energies of all the point 
charges, which originates in the 6(r) term in Eq. (38). Sub- <, 
tracting this off will give a more revealing quantity. The ^ 
self energy of one charge is proportional to E| V g = (8ue /Z V) 
/ dk/k, and the self field-energy for 2N of them is an additive 
quantity, so 


(f) 


self 


16 rrn e c 
o 


r 


dk 

k 


and 


<® 2 > - <*\elf " 


2 e 

*4 


r c 


JO 


dk 


k(l + k 2 ^) 


(45) 


Note that Eq. (45) has now a small k divergence, which must be 
associated with targe spatial separations* This has no analogue 
in the three-dimensional theory, even though (23) is essentially 
the same as the three-dimensional expression. This pathology 
has its origin in the fact that the Coulomb potential (15) has 
a divergence at infinity in two dimensions, but the three- 
dimensional expression does not. The divergence strongly sug- 
gests that we may anticipate other unexpected small-t diver- 
gences , and that at times we may want to limit the discussion 
to a large but finite plasma volume V = L ? , to provide an effec- 
tive cutoff in k space at k " 2 tt/L. 


An interesting ratio is the ratio of the electric field 
energy density to the thermal energy density 2n Q 0/£. Dividing 
(44) by 167 Tn o 0 /& and cutting the integral off at x min , * *d € l/ " 



MAGNETIZED PLASMA MODELS 


449 


<!!> 

I6nn o 0/^ 


8tt n o ^ 


X*/ 


( V 


x . ) 

min. ' 


--i-eue 1 /*, i 

This is not small small compared to unity for any However, 
if we compute the comparable ratio for the part of <E 2 > that is 
to be associated with interaction energy. 


<3 2 > - <g 2 > 

I6nn 0 9/r £ " 


(46) 


As long as the right-hand side of (46) is small compared to 
unity it makes sense to speak about the field fluctuation energy 
as being small relative to the kinetic energy, which is the 
essential approximation that has been made. But it is clear 
that such formulas as (46) lead to grave difficulties if one 
tries to pass to the limit L Whether these are difficul- 
ties with the perturbation expansion (17) , or correspond to 

more basic questions connected with the existence of the thermo- 
dynamic limit in two dimensions is something that at the time 
of this writing, apparently nobody knows. Definitive answers 
to these questions will Involve investigations of a consid- 
erably higher degree of rigor than the one presented here. In 

the meantime, one wants to be especially cautious of attempts 

to evaluate extensive quantities, sujjh as the total energy, by 
integrating intensive ones such as <E 2 >. The circumstance may 
arise (it does in one dimension [15]) in which small terms of 
G(l/V) may be misled in the perturbation theory, which give no 
contribution to <E 2 > in the limit, but which may nonetheless 
give finite or even divergent contributions to integrals over 
the entire volume. 

(H) Equation of State , Collapse at Lew Temperatures 

As a next comment on the equilibrium theory, let us remark 
on one rather remarkable property of this system, namely that 
its equation of state is exactly calculable. This has apparently 
been discovered at least twice [16,17]. It is not a complicated 
manipulation [2] to reduce Eq. (4) to 



450 MONTGOMERY 



In (48) , / d^ N ^ X means the integral over all phase space co- 
ordinates except for x. and x . . The is over all pairs of 
both signs. 3 1<J 

Equation (47) applies to any force law <\>±j . The unique 
feature of the two-dimensional Coulomb potential is that 

(xj_j) - -2e i ej/£ > independently of x-^ and Xj . What re- 
mains is just the normalization integral for the two-body dis- 
tributions, and the problem has been reduced to counting par- 
ticle pairs. This gives, after a little algebra 



(49) 


The equation of state, while becoming that of an ideal gas for 
high temperatures (e « 1) , has the feature of predicting a 
collapse (p < 0) at low temperatures, corresponding to e - 4. 
This is nothing remarkable, since it has its origin in the 
fact that the ion-electron interaction (and thus the total 
energy) is unbounded from below. In fact, one sees from (3) 
that when x-^ is close to Xj , Z contains a factor 



1 2e e /A0 

~ x. . dx. . x. . J 
J ij aj 


If and ej are of unlike sign, this will diverge when -(2 2 /£9) 
+ 1 < -I, or e > 2, If more than two particles are considered r 
the divergence may even occur at higher 0 (see Knorr [18]). 

Some lower bound on the interaction would be required to elim- 
inate this collapse, as quantum mechanics is required in three 
dimensions . 
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(I) Summary 

Starting from the Gibbs distribution for an equilibrium 
system of electrostatically interacting charged rods aligned 
parallel, approximate expressions for the phase space proba- 
bility distributions have been derived by an expansion in the 
plasma parameter (number of rods per Debye square)* The pro- 
gram has been carried out in parallel to the familiar [10,13] 
expansion in three dimensions* These approximate probability 
distributions permit the calculation of ensemble averages of 
stationary fluctuation quantities and their spectral densi- 
ties; these are independent of external magnetic field strength. 
The thermodynamic equation of state has been noted to be exact. 
Anomalies which may possibly result from taking the "thermo- 
dynamic" limit (N 00 V TO ) have been noted. Up to this point 
the theory is the same with or without an external magnetic 
field. 


II. THE TWO-DIMENSIONAL GUIDING- CENTER PLASMA 
(A) Some General Remarks on Transport Processes 

A central concern in statistical physics has always been 
the matter of transport properties [19]. That is, gradients 
in the macroscopic hydrodynamic parameters will result in the 
net transport of mass, momentum, energy, electric current, etc., 
across an imaginary surface in a gas, liquid, or plasma. This 
transport comes about as a direct consequence of the microscopic 
interactions among the particles, and is inevitably a problem 
in non- equilibrium statistical mechanics. 

This problem has been of acute concern for strongly mag- 
netized plasmas, because many controlled fusion schemes pro- 
pose to confine a very hot plasma away from material walls by 
means of externally- imposed dc magnetic fields. Often, actual 
machines appear to confine plasmas less well than they are sup- 
posed to, and this "anomalous" transport of plasma has been a 
matter of concern for many years. It was addressed early by 
Bohm [20], disappeared from sight for several years, re-appeared 
about 1960 [21,22] when controlled fusion research became 
public knowledge, and has been the motivation for many of the 
exotic instability and "weak turbulence" calculations which 
were a major pre-occupation of plasma theory in the 1960 t s. 

The subject is far from exhausted, though it has become dif- 
ficult of access by virtue of being surrounded with hundreds 
of rather diffuse calculations. 
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Actually, the transport may not be so "anomalous 11 as it 
might appear at first sight* Most of the more sober calcula- 
tions [23] .have started from a kinetic equation . By kinetic 
equation, one means a differentio- integral equation such as 
Boltzmann 1 s equation, which will advance the time- dependent 
analogues of f lj0 and f 0jl of Eqs. (9) in time* A crucial 
ingredient of any kinetic equation is a collision term in 
which the irreversibility originates, which is a functional 
of the one-body distribution (f x 0 or f 0 x ) , and into the 
derivation of which goes some mi£roscopi£ model of the par- 
ticle-particle interaction process* Different forms of the 
collision term arise from different microscopic models [24], 
Thus uncorrelated, strong, two-body collisions lead to the 
Boltzmann equation, and weak two-body collisions lead to the 
Fokker-Planck equation. In the plasma case, however, most of 
the microscopic models used to describe the interactions of 
charged particles have not included the effects of an exter- 
nal magnetic field on the particle-particle interaction pro- 
cess. Putting in an external magnetic field renders the 
strongly- interacting two-body problem insoluble, so there is 
at present no analogue of Boltzmann f s equation for the case in 
which an intense, external magnetic field is present. In a 
weak-interactlon model of a plasma, a kinetic equation results 
[25-27] but is so complicated that even the derivation of it 
is quite strenuous, and heroic efforts will be required to 
extract many predictions about transport properties from it. 

From the foregoing remarks we may conclude that the sub- 
ject of transport properties in a strongly-magnetized plasma 
is still very much an open one* One expects significant 
gains in understanding only to the extent that clean and bold 
simplifying approximations can be found* A recent breakthrough 
in this direction was due to Taylor and McNamara [28] . They 
introduced a two-dimensional model of a plasma consisting of 
charged rods aligned parallel which are free to move perpenj- 
dicularly to a strong dc magnetic field. They further simpli- 
fied the dynamics by introducing the guiding- center approxima 
tion for the particle motion. Discussion of the various prop- 
erties of this model will be the subject-matter of this chapter. 
Various features of the model can be generalized to three di- 
mensions, but discussion of these is better deferred until 
after the two-dimensional model has been presented* 
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(B) The Guiding- Center Plasma ("Zero Gy to- Radius " Plasma) 

The physical system is the same as that described in 
Section I: a collection of N very long rods of charge +e, 
length &, mass m-£, and N more of charge -e» length mass m e 
aligned parallel to the z-axis of coordinates, and interacting 
through the Coulomb potential, with an interaction potential 
given by <f>„ m = In Ix^ - .xj for any two charges. 

The equation of motion for the velocity of a given charged 
rod of charge-to-mass ratio e j/ m j 

„ !i (l + | xt) (50) 

dt> V c J 

jfhere v* , E are two-dimensional vectors lying in the xy plane, 

B = Be z , with B - a constant, and E is evaluated at the loca- 
tion of the particle. E 55 -V0 obeys the two-dimensional Poisson 
equation. 



(51) 


where the sum runs over all charges of both signs. Equations 
(50) and (51), together with dxj/dt = and the appropriate 
boundary conditions, constitute a complete dynamical descrip- 
tion of the system. This dynamical description is still. quite 
complicated, and we shall have more to say about it later. At 
present, we wish to discuss an even more simplified version 
of the dynamics, namely, the guiding -center approximation * 


The reader is referred to Northrop [29] for a discussion 
of the guiding center approximation; in practice, it amounts 
£o neglect ing^the left-hand side of Eq. (50), so that 
E + (vj/c) x B ^ 0, or in this two-dimensional geometry, 

v. cs c l?(x.,t) x 3/B 2 (52) 

J 3 

The Vj in (52) is to be distinguished from the vj in (50) ; 
(52) describes the motion of the guiding center of the par- 
ticle. This guiding center motion has superimposed on it the 
"fast" gyration of the particle with characteristic gyro- frequency 
Q = e B/m^c and characteristic radius of the order of - 

i -i 3 3 

/O/m V/ftj , for a particle of thermal energy 9. The "fast" 
motion is averaged out in writing down Eq. (52). The conditions 
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for the validity of Eq. (52) as a satisfactory approximation 
are that, most importantly , the characteristic time T and 
length A over which E(x,t) varies significantly must satisfy 
the inequalities 


X » r 

A 

n j T >> i 


(53) 


The actual microscopic electric field in the plasma will be 
de- comp os able into components which do and do not satisfy the 
inequalities (53) . Only the former will be treated accurately 
by (52)* At the present time, no satisfactory quantitative 
estimate exists as to the errors which result from this ap^ 
proximation. v 

(C) Connection with Ideal Hydrodynamics 

The replacement of the differential equation (50) by the 
algebraic equation (52) is a considerable simplification. It 
also has the unexpected benefit of leaving us with a Hamiltonian 
system, the Hamiltonian of which is essentially the same as 
the potential energy of the original system! To see this, we 
digress briefly on the hydrodynamics of incompressible, inviscid 
fluids in two dimensions. In the theory of parallel "line" 
vortex motion [30], it is shown that the fluid velocity field 
which results from a collection of line vortices of vortex 
strength Kj = Kje z located at positions xj is 



(54) 


The vortices move along streamlines, so in particular the 
velocity of the ith vortex is 


£ = > JL x JJ. 

v i 2 tt * 2 


Now consider (50) , (51) , and (52) . The electric field seen 
by charge i is 
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so its velocity is 


(56) 


^-1 

X. . 

s' 

“Z, i b 

2 X 

B 


X id 



(57) 


If we identify ( 2 tt) — 1 times the vortex strength, K,/2 tt, with 
-2cej/£B, then the two mathematical descriptions become iden- 
tical . We may use language appropriate to the guiding center 
plasma interchangably with language appropriate to line vorte> 
motion, providing that 


K. = -4ttc e / & B (58) 

d u 

That the system is Hamiltonian with Hamiltonian 


H - -h E K i k j l 

i<J 

can readily be seen by picking canonical coordinates 

[ *j = (x j’ y j )] 


for each vortex (or charge) . 
is the same as 


P 


i 



i 


% 


sgn K.) 

<■ 

It is readily verified that 



(60) 

(55) 


(61) 


Thus the canonical coordinates and momenta are essentially 
just the Cartesian spatial coordinates x^y^ of the charges 
(or vortices), with a change in sign of for the negative 
vortices* 

Since H is not explicitly time dependent, it is a constant 
of the motion. The canonical distribution is again given by 



456 


MONTGOMERY 


Eq. (±) , since H is up to a constant factor, just the inter- 
action energy ^ of all the pairs of charges. All the 

stationary thermal equilibrium correlation functions and spec- 
tral densities calculated in Section I are identical to the 
corresponding quantities for the guiding center plasma* Of 
course time dependent quantities will evolve much differently, 
and we pass now to a consideration of some of these. 

(D) Formulation of the Problem of Spatial Diffusion in the 
Green- Kubo Picture 


A measure of the rate of spatial transport can be obtained 
in the following simple situations. Suppose a very tenuous 
distribution of labeled particles, so rarified that their inter 
action with each other is negligible, is released near x - 01 
As time goes on, due to the fluctuating microscopic electric 
field, the charges (they will hereafter be called "test 1 ' par- 
ticles) will be moved away from the origin. The position after 
time t for one of them will be 


$ = f v( T ) a t 

'0 


where v(t) is the velocity. By (52) 3 this is 



(62) 


(63) 


where E(t) = E[x(t) ,t] is the electric field evaluated at the 
Instantaneous test particle position. 


The path (63) is erratic and complicated, and to evaluate 
it would involve knowing, among other things, the exact micro- 
scopic electric field E(x,t). An ensemble average prediction 
is the best we can hope for. So we imagine a set of similar 
plasmas, distributed according to the Gibbs distribution (1), 
and release a collection of test particles near J = 0 at t * 0 
in each one. Hereafter we address ourselves to the question 
of calculating ensemble averages over this ensemble [31]. 

. 

Since <E> = 0, by Eq. (33), we have at once that <x> 53 0, 
which simply says that particle transport is isotropic. A 
measure of the amount of diffusion is <% 1 > f 0; thus 
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Auto -correlations like <E(Tj) * E(t 2 )> play a central role in 
what follows ♦ If the ensemble is stationary, they will depend 
upon t 2 - f i only, so that <E(T X ) • E(t 2 )> = <E(0) • E(t 2 “ T i^ >J 
and it is expected that for all physical systems, <E(0^ E(t^> 
will -* 0 as T ■+ ± 00 • It is also to be expected that <E(0) E(t)> 
is an even function of T. We shall assume so. 

Defining X 2 - T, = x, (t 2 + T t )/2 5 T, the integral (64) 
can be rewritten as 

(x 2 ) = ~ J*dT J* <f(0) * *(T)> (65) 

The region of integration in (65) is as shown below (Fig. II- 1) . 



Figure II-l. Regions of integration for Eqs. (64) and (65). 

For large t, the non-vanishing part of the 
integrand will lie in the neighborhood of 
the T=0 axis. 


Now if <E(0) * E(T)> goes to zero sufficiently rapidly 
(more rapidly than 1 /t) for larg£ T, we^may consider the case 
when t is >> the time in which <E(0) • E(t)> goes to zero, and 
write, from (65) , 
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#->2\ 

<X ) 



<3(o) • 3 (t)> 




( 66 ) 


The ratio <3 2 >/2t thus becomes time- independent for long times 
and this ratio D, 

2 r 00 

0 - / <3(0) • f(T)> dT (67) 

B Jo 

is called the diffusion coefficient . It is generally thought 
of as a good measure of the diffusiveness of a medium. (It 
is also frequently defined by <5 2 >/2td, where d is the dimen- 
sionality of the system) . 


Recently a lot of interest has arisen in the possibility 
that <x 2 > may grow faster than t as t *» in two ^.imensigns 
for a variety of systems [32]. For example, if <E(0) * E(t)> 
were to fall off as 0(1 /t) for large x, detailed examination 
of (65) reveals that <t 2 > - t In t for large t. This has been 
hard for many people to accept, but it has also shown up now 
in the plasma case. The reason the <x 2 > ^ t behavior is so 
well accepted is that it can be derived regardless of dimen- 
sionality [31] , from a random-walk model in which a random 
walker takes a large number n of uncorrelated steps per unit 
time, each step of very small length a, with n<a 2 > = const.. 

(D is in fact, just proportional to n<a 2 >.) This is true for 
virtually any distribution law for the steps. If <i 2 > turns 
out to grow faster than t for systems of dimensionality less 
than three, then it means that the venerable random walk model 
is probably less accurate than we had supposed for these sys- 
tems . 


In any case, we will be concerned in the following pages 
with attempts to calculate the function 

2 f t 

D(t) s ~ <f(0) • 3 <t)> dT (68) 

Jo 
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for the guiding- center plasma whose dynamics are given by (51) 
and (52), If the integral converges, it will be called the 
,r dif fusion coefficient' 1 , 

D 55 D(°°) (if it exists!) . (69) 

But if the integral D(t) fails to approach a limit, it still 
is an interesting measure of the way the test particles diffuse. 


(t 


(E) Evaluation of 


<E( 0) • E(t)> d'l and D 


o 


It is convenient t£ evaluate (68) in a large but finite 
volume V = L 2 . E(t) - E[x(t),t], and the exact electric field 
due to the 2N charges is 


2(3,t) = - £ 


k 


4 nil? „ 

• — p~ P.W e 

k k 


where the Fourier transform of the charge density is 


(70) 


1 j-* -ik’x > 

p^(t) B — dx e p(x,t) 


(71) 


Using (34) , and inserting the result into (70) gives 


a-* r -1 r -1 e « 4 nil? [x-3 (t)3 

S(x,t) = - y L *f " 


2 . 

XV 


(72) 


* A 


where the §um is over all charges of both signs, and the 
values of k have components which are integer multiples of 
2 tt/L (we impose periodic boundary conditions). 


E(t) = E[x(t), t] is then 

. r-, 4ne it iit-[x( T )-x.(T)] 

i(r) = V - - g - a — e 0 

fr 1 ik j ev 

k,J 


/ 


( 73 ) 
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where x(t) is the trajectory of a particle which starts at the 
origin: x(0) = 0. For shorthand, (73) will also be written 

as 



in obvious notation. 

The problem is then one of evaluating 



This result is exact, and we cannot go further without 
an approximation or an hypothesis. One which can be made [28] 
but which some day will have to be given careful scrutiny is 
to neglect the correlations between the positions of the test 
particles and those of the background plasma, the (t) ♦ That 
assumption reduces (75) to 


<3(0) • 2(t)> 5 




Further, noting that^E* (t^ = 0 unless k^ + k 2 = 0 

for a spatially-uniform ensemble, and that E £(t) = E^(t) since 

->• 

E is real, we have, from (76) , 


<2(0) • 1(f)) = £ («>> • y T j)^e ik ' x(T) ^) (77) 
2 

The next step necessarily involves expressing the two 
terras on tlje righ^-hand side of (77) in an approximate way in 
terms of <E(0) • E(r)>, to provide an equation to be solved 
for this quantity. The last term can be written as 


<exp[-ik • x(t)]> 
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where E(m) is the electric field between times (m — 1)t/M and 
mT/M* At is t/M. 


The computation of the expectation value in (78) involves 
not just the probability distribution of the electric field 
seen by the particle, but the joint probability distribution of 
£he electric field at a whole sequence of times, E(l) , E(2) , 
E(3). . , If we ca^l this joint probability distribution 
,[E(1) ,E(2) , . • . E(M) ] , then formally 

<exp[-iS • x(t)]> 


= lim 

M*+co 


r 

d3(l) ... d2(M) P[2(l) ... 2(M)] 


exp 


icE x 


3 


M 


fit 


B 


Y 2(m) 


m=l 


(79) 


It has not been possible to establish an expression for 
P[E(1)- ... E(M) ] from first principles [33] . A reasonable 
conjecture is that it is a jointly normal probability distri- 
bution [34,33]; that is, 


P[2(l) ... 2(M)] = T| exp 


M 

•E 

L i,J=i 


n 


a..:2(i) E(j) | (80) 


where n is a normalizing constant and the are constant 
2 x 2 positive-definite dyadics. A generalization of the cen- 
tral limit theorem, for example, makes (80) look extremely 
plausible [35]. In tjje present context, the theorem s^ays t£e 
following. Consider E (1) >E( 2) , • • « to be the sum Zj [E^ (1) ,Ej ( 2) , • * ■ ] • 
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Here Ej (n) is the electric field at the test particle produced 
by particle j in time interval n, and the Zj runs over all the 
plasma charges. Then, under very weak conditions oij the moments 
of the probability distributions of^the individual Ej (n) , and 
the additional assumption that the Ej (n) are uncorrelated with 
each other for different j , the probability distribution of the 
Zj (where j goes from 1 to 2N) converges to (80) for N + 

This neglect of correlation between the trajectories of the 
plasma particles is the key approximation, and is made quite 
plausible by the weakness of the pair correlation [Eqs. (21) 
to (25)] for small plasma parameter. 


It is shown in Appendix I that the result of substituting 
(80) into (79) is 

<exp(-iic • x(t) 3 > 


lim exp 
M-4os 


M M 


■t[ E <%) %». y tt H 


m l =1 m 2 =1 


= exp 


2 ,_2 


c k 


2 B 


T T 

J dT 1 J dT 2 <%) E(t 2 )>: (k X B)(k x B) 


exp 


77 J dT l / dT 2 < g < T l> • ®<V> 


0 o 


(81) 


Inserting (81) Into (77) yields the result that 


< 2 ( 0 ) • 2 (t)> 

Z (V <0) ■ V T> ) exp ' 


4 B 2 






Jo Jo 


d T 2 <&(')) • 3( t 2 )> 


( 82 ) 
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Equation^ 82) a satisfactory starting point for cal- 
culation of <E^0) • E(j)> as soon as we have an acceptable ex- 
pression for <E^(0) • E+(t) . First, note that 


<C<°> • ^j t )> 

k k 


E 


16 

2 2 2 

irf 


ig-[t(0)-x 4 (x)] 


> 


(83) 


V 16 "S % 

f; iW 



r T 8 (s) x 8 
x a (°) - 8/o) - « « 

Jo — 


By E. (t) , we mean the electric field seen by the Jfcth particle, 
of charge e^. In (83) , we can treat the £th particle in the 
same approximations we treated the test particles in Eqs . (76) 
through (82) , getting 

<C<°> • 

k it 


_ 16 n e . e . x 

^ — P ~ ^ - — /exp ik • [x^(0) - *,(0)]^ , 


A,t 


2 2 
i k V 


/ let X 3 f 
/exp ^ • j 

\ B J o 


e a (5) *5 


’ <V (0) • V 0) ) e:ip {■ 71F I dTl I 4 " 2 <S(,1> ' ?(T2,> } 


( 84 ) 
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Since the E x B drift is independent of both _£he mass and sign 
of the charge, the statistical properties of E^(£) will be the 
same for all 

Inserting (84) into (82) gives 


<£(o) • t)> 


-E (IS |2> ex,). J | *2 (S( V %») 

£ 1 o ° J 

Since <|eJT| 2 > is a stationary quantity obtainable from 

the methods of Section I, Eq * (j£5) is ijow a closed equation 
that can be used to determine <E(0) • E(t)>. Defining 
R(t) = ( c 2 /2B 2 ) • 

f V f T dT q(T 2 - T ), where q(t) = — <f(0) * ]£(t)> 

J O X J o 

dR(T)/d T = (c 2 /b 2 ) f T Q(5) dg, and d 2 R/dT 2 = c 2 Q(t)/b 2 , which 

Jo 

means that (85) can be expressed as 


2 2 

2 % — £ = Y<\K\*> exp [-2k 2 R( T )] (86) 

c dT^ K 

i 

From (65), R(t) has the physical interpretation of being 
<£ 2 >/4* where <x 2 > is the mean square spread after time t of 
the distribution of test particles. 

Equation (86) has a first integral 

E<i 

k 





MAGNETIZED PLASMA MODELS 


465 


where the fact that dR(0)/dT = 0 has been used to evaluate the 
constant of integration. To proceed further with Eq. (87) re- 
quires numerical integration, if R(t) is desired. This inte- 
gration is discussed in a later section. We can extract the 
long time behavior analytically, however, by noting that 
R(t) + °°, so that 



approaches the limit 


D(«) = 2 



< ^ 1 2 > 



Up to a factor of /2, this is the result of Taylor and 
McNamara [ 28] . 


( 88 ) 


(89) 


Several remarks are in order concerning (89) before we 
pass to the details of its evaluation. First, it varies as 
1/B, confirming the elusive conjecture of Bohm [20K Secondly, 
it will contain a factor /e, since it contains £ <]E^| 2 >/k 2 , 

k 

and will vanish in the Vlasov limit (z 0). Finally, the ex- 
pression will diverge in the infinite volume limit, since the 
minimum value of jk] is 2^/1, 9 and a factor 1/k 2 occurs in the 

This is a reflection of the fact that the long-wavelength, 

small -k contributions to (89) are the dominant terms, in con- 
trast to "classical' 1 diffusion models based on random two-body 
encounters. 


We may estimate (89) for a larg^ volume V - L 2 by replac- 
ing the sum over discrete values of k by an integral over con- 
tinuous k according to the prescription 
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Trace 


Sg(fi) 


= 2 TT 


k dk 


JB. 

Itx 


i * k s 


(90) 


where use has been made of (43) . The integral in (90) is to 
be cut off at k min = 2tt/L from below. The result of using 
(90) in (89) for the diffusion coefficient is 



r 

r 00 1 

2c 


dk 

B 

L‘ J 

k . k(1 + k2 
min 


2c 

B 



Iaj (L/2 rr 



(91) 


[jjrtt. (L/2 tt \ d )] % 


The Bohm result, never justified in detail, was D(°°) 


JL 

16 eB* 


Notice, however, the following limitation in the deriva- 
tion of Eqs. (89) and (91). It was arrived at by assuming 
that the terms involving exp(-2k 2 R) in (87) had become « 1* 
Since “ 2TT/L and R = <x 2 >/4, this is equivalent to 


2 




» 


1 


or <5 2 > » L 2 /2 tt 2 H L 2 /20. However, when <x 2 > becomes i L 2 , 
the test particles will have encountered the boundary once, 
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and in view of the periodic boundary conditions, the diffusion 
model with its monotonically increasing <£ 2 > will have broken 
down* Thus there is a relatively narrow temporal interval, 

T 2 2 

« 2b D(®) « IT (92) 

in which (89) and (91) may be expected. to apply* Both before 
and after this interval, the development of <x 2 > will be more 
complicated. 

Finally, we remark that if one attempts to apply (86) to 
unboimde d tw o-dimensional plasmas, a variation [36] of 
<x 2 > - t/lrTt occurs for large values of t. 

(F) Liouville Equation* BBGKY Hierarchy* and the 
Time-Dependent " Vlasov " Limit [12] 

The diffusion calculation could proceed in a relatively 
unsystematic way, with the assumption (80) as the key ingredient 
For other purposes, a more systematic formalism for the kinetic 
theory of the two-dimensional guiding center plasma is desirable. 
Proceeding in parallel to the well-known program of the non- 
equilibrium BBGKY [2-6,11] hierarchy, it is useful to construct 
a description based on the Liouville equation . 

The phase- space probability distribution of any Hamiltonian 
system obeys a Liouville equation, which simply expresses the 
constancy of that probability distribution as seen by an ob- 
server who moves along any phase-space trajectory traced out 
by the system as it moves through its phase space. Since we 
have seen the system to be Hamiltonian in Eqs. (59) to (61), 
we can write down the Liouville equation immediately as 



where D « D(xj ,x 2 , . * * , t) is the probability distribution of 
all £he particles in their phase space. Noting that dx i (t)/dt 
= c E(xi,t) x B/B 2 , (93) becomes 



( 94 ) 
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where tlCx-^t) = + ^ ^ ^ ^ ^ 1 X i " X j I > and t ^ ie 

9 X i j?fi 

summations run over all charges of both signs. The single 
most important solution of (94) is. of course, 

D eq. = ex P ( “E <P itj /9)/I dX exp ( - £ tp^/e), though for a 

non- equilibrium situation, D will not be D g . We shall gen- 
erally assume, however, that the various correlation functions 
have the same orders of magrvitude as the thermal equilibrium 
correlations of Section I, We also assume D to be symmetric 
under interchange of like particle coordinates. The reduced 
probability distributions are defined in analogy to Eqs. (9) 
as 


f „ = V 
S.jS A 

i e 


s.+s 
i e 


^v 1 


where it is to be noted that 



+ l 

x 

s. 

1 



is a function of position space coordinates only (the dimen- 
sionality of the phase space is now 4N) . 


Integrating (94) over a big box of volume V = L 2 gives: 
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(96) 

This is the non- equilibrium BBGKY hierarchy for the system, 

and is the analogue of (13)* By E^ Si>Se \ we mean - — 

J 3xf 

J 

^ s i ,s e^* w ^ ere $ (spS.) is the £ (f> for the first Sj_ ions 

3 i<j 3 


and the first electrons, 
e 


Similarly E;£ s i ,s e> 


-3 


Equations (94) and (96) are not as complicated as they 
look, and we can profitably introduce abbreviations for the 
linear differential and integral operators which appear in them. 

Thus, (94) will sometimes be written as ( I HR )D - 0, and (96) 

' ot o / 



H( V 8 e> 



f r ii T ie \ 
f L + L i 
l s« s , s | 

V i e i e/ 


's.+l>s 

i > e 


+ 



+ L 


ee 

3. S 

i e 



v 1 


in an obvious notation. H 0 is the 


as 
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Liouville operator for all 2N particles, and H(s ±9 s e ) is the 
Liouville operator for si ions and s e electrons. 


The BBGKY hierarchy (96) is analogous to the well-known 
hierarchy for ordinary plasmas but it has important differences 
which are worth mentioning. It is interesting to consider the 
"Vlasov limit" in which the correlations vanish; the corres- 
ponding thermal equilibrium limit (see Section I) is that in 
which e -*• 0. Thus we write 


V S e 


6 . 

1 


n f i o^ t} 

j=i 1,v 3 


s 

n f 

k=l u ’ 3 


and all of Eqs. (96) 


become satisfied if 


9 f l,0 ,-*i . . c 1 ^ 5 f i.0 (x i ,t) 

at K i’ > b 2 •* 


= 0 


a f n , „ c 


f Cx 

o,r l- 


3 x 1 


(97) 


,t) 


= 0 


where 



Equations (97) and (98) are the analogues of the classical 
Vlasov-Poisson system for the two-dimensional guiding- center 
plasma. The spatially-unif orm state, f ly0 = f 0>1 - 1 and E = 0, 
is obviously a solution. However, the analogue of the Landau 
problem, in which one linearizes about the spatially uniform 
state, is apparently insoluble; for the first non-vanishing 
terms in (97) are intrinsically of second order in powers of 
the departure from spatial uniformity.- 

Equations (97) , (98) are completely equivalent to the 
equations of motion for incompressible, inviscid flow. They 
can be made to look simpler by subtracting Eqs. (97) from each 
other to get 
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St 


+ 



where 


dx 


0 


■— • E - 4 TT p 

3 * 


(99) 


( 100 ) 


Very little information exists on the solutions of the 
system (99) , (100) , although a number of numerical simulations 
exist which are relevant [37-39] . Taylor and Thompson [40] have 
argued for the existence of a large class of oscillatory so- 
lutions of thi^ system with a frequency which vanishes as the 
amplitude of [e| 2 vanishes. This derivation, however, is not 
wholly transparent. Any oscillations which do exist will be 
intrinsically nonlinear. 


Another peculiar feature of the hierarchy (96) is the 
apparent absence of any tendency toward thermal equilibrium. 

For the conventional hierarchy for a normal plasma, assuming 
an initially spatially uniform state leads to rapidly-developing 
pair correlations. These in turn relax to functionals of fi o 
and f 0>1 which, when substituted into the equations for df x ’/Bt 
and Bf 0 l /3t, lead to the Balescu-Lenard equation [7,41-43].* 

From this, an approach to thermal equilibrium can be proved. 

The situation is quite different for the guiding-center plasma. 

It is clear by inspection that the completely uncorrelated state, 

i S e 

^ /HS 


f l,0 ( V f 0,l ( *k )With fl ’° 


■ 0,1 


in which f = R II 

S i ,S e j=l k=l 
is a tijne-iTidependent solution of (96) I It is not known whether 
there is any sense in which the solutions of (96) relax to the 
thermal equilibrium values (11) . This is one of the really 
mysterious parts of the theory. A number of time independent 
solutions with non- thermal equilibrium correlations for (96) 
can be constructed, but there is no particular reason for pre- 
ferring one over the other. 


- !• 


(G) Calculation of the Electrical Conductivity by the Kubo 
Method [44,45] 


A calculation of a second transport coefficient, the elec- 
trical conductivity, proceeds most naturally from the Liouville 
equation, Eq. (94) will be abbreviated as 



H 


D 


0 


(94a) 
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The addition of a spatially uniform electric field E Q e 
modifies (94a) to read 



where the right-hand side Eq. (101) defines the linear op 
erator H . We now regard E Q as weak , and seek a perturbation 
solution 1 Co (101) in the form D = d(0) + d(D + .... where 
the superscript indicates the order in . Thus ' obeys 



d<°) 


s= 0 


and D 


CD 


obeys 


— + D^ 1 ) = - exp(-ia)t) H 

o 1 


( 102 ) 


(103) 


The appropriate solution to (102) is the Gibbs distribution 


= T| exp (- £ (p^/G) ^ 10A) 

and the £ cj^j runs over all pairs of particles of both signs. 
i<j 

As usual, n is chosen so that J D^) dX = 1- Since (102) is 
simply a statement of the constancy of d( 0) along a trajectory 
in the total phase space, it is just a statement that the 
formal solution 

D (0) (t) = e tH ° D { °)(0) 


(105) 
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is just a function of the constants of the motion [6]* £ 

i< j 

has already been noted to be such a constant, and (104) is the 
solution appropriate to thermal equilibrium. 

-y 

We imagine the external uniform field to be. turned on 
abruptly at t = 0. The solution to (103) can be found by in- 
tegrating the equivalent equation 

± (.“« >><») . - % D<°> (106) 

between t - 0 and t = t. Setting 8 0 at t ■ 0, noting 

that e“ tH o = d( 0) , and rearranging the integral gives 


D h) „ _ 


d T e i<uT e ° 


(107) 


Now for any dynamical variable whose value is determined by 
the phase space coordinates (A, say) , the expectation value 
<A> = J A D dX. For variables which have expectation values 
which vanish in thermal equilibrium, we have the result that 
to first order in E 0> <A> = J A d( 1) dX. Introducing the 
explicit form of H 1# we "have 



1 


We now note that the effect of the operator e ° applied to 
any function of the coordinates obeying Liouville’s equation 
is to trace those coordinates back t units along the phase 
space trajectory and assign them their values there. Thus, 

“ TH ° -y-J -yi 

for instance e • v. = v^(-t), where v.(-t) is the velocity 

j j J 

that the jth ion would have had x seconds- ago if it now has 
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velocity v* = (0) . Therefore 


exp 


(-tH d ) E 1 D (0) = - } ^ e pJj(-T) - v*(-t)J D (0) 


(109) 


,( 0 ) 


since D is constant along the phase space trajectory. In- 
serting (109) into (107) and the result into the expression 
for <A> gives 


<A> 


= g*p( 


-iuot) ^ 


dT exp (iicx) 


(110a) 


<A T e • [v*(-t) - v®(-T)]> 
(} 


We have not yet committed ourselves to an explicit expres- 
sion for the dynamical variable A. To calculate the electric 
current density produced in response to the electric field E 0 , 
we pick A a function which makes <A> equal to the volume cur- 
rent density <j>. We must remember to apply this procedure 
only far away from the plasma boundary. Thus 

A •+ ^ = X e (v* - v\)/ II? 

5 J J 

Noting that <v(0) * v(t)> will be an even function of t for 
all particles of both signs, 

<2> = 

U) 


e ex-p(-iqjt) 


XL 2 e 


o 

-L 


dT e 


iu)T 
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<[Vj(0) - Vj(O)] [v^(T) - v^(t)]> . (110b) 

i?& .+ 

The coefficient of E Q exp (-icot) in (110b) can be identified 
as the conductivity tensor 8(b)) : 


a(«) 


s f iujT 

2 dT e 

£9 L J o 

^ <(v^(0) - Vj(0)) [v*( t) - v®(t)]> 

t A 


(111) 


The terms in (111) with j ^ % and not of the same species can 
be shown to give a contribution which is down by a factor of 
0(1/N) from those with j = £. Therefore 5(a)) becomes, since 
all the j * l terms of the same species are equal. 


a(u)) = 


2Be 


e iujT <v(0) v(t)> 


2 n e 2 c 2 
o 


x e b 



( 112 ) 

dT e iUJT <3(o) f( T )> x i 


The <E(0) E(t)> is the same quantity whose trace appears in 
Eq« (64), and the same methods can be used to calculate it*_^ 
Under the apparently straightforward assumption that <E(0) E(t)> 
is diagonal (for it to be otherwise would imply a preferred 
direction in space) , we have 


ct((u) = S cr(iu) 

where 


o(u0 = 


o 2 2 

2 n e c 
o 

IBB 2 


Q(t) e d<WT dr 


(113) 
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with Q(t) given by j <E(0) • E(x)>, as before. 
In the dc limit, to = 0, 


a(0) o 


2 2 

2 n e c 
o 

JL 9 B 2 


f 


Q(T) dT 


n e 

0 

1 e 


D(«) 


(114) 


showing a direct proportionality with the diffusion coefficient 
D(°°) of Eq. (89). Such a proportionality between the conduc- 
tivity and diffusion coefficient was called by Kubo [44] a 
"generalized Einstein relation", for systems whose dynamics 
are governed by Newton’s laws. But it was far from obvious 
that such a proportionality would exist for the guiding-center 
plasma. 

The ac conductivity requires a Fourier transform of Q(x) , 
which is not obtainable analytically in closed form. For that, 
one needs to solve Eq. (87) numerically. This has been done, 
and the details are reported by Montgomery and Tappert [45]. 
Typical results are shown, in dimensionless units whose inter- 
pretation we can ignore here, in Figs. II-2 and II-3. The 
same remarks about divergences for L “► 00 apply as were ap- 
plicable in the discussion of the diffusion coefficient. 

(Note added in proof : Recent numerical simulations by 

G. Joyce and the author have indicated that the conductive 
behavior of the guiding center plasma may be somewhat more 
complicated than the present development would indicate, A 
transverse electric field induces large scale vortex motion of 
a non-local character in a wide variety of simulations. We are 
also indebted to Dr. Leaf Turner for pointing out to us certain 
shortcomings of the present derivation of a local conductivity.) 
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Figure II-2. Numerical solution of' Eq, (87), from Mont- 
gomery and Tappert [45], Dimensionless units 

2 L 

are: k/Kp, R/2Kj), t[8TTnoce£ ]/£E. Notice 
that R increases linearly with t almost from 
the beginning. 



FREQUENCY w 


Figure II-3. • Real and imaginary parts of the conductivity 
from Eq. (113). The units of a are noeceV^B 
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(E) The " Negative Temperature” Instability for the 

Guiding- Center Plasma 

A rather unusual feature of the guiding- center plasma is 
that for interaction energies greater than a certain 

amount, statistical mechanics predicts no stable, thermal equi- 
librium, spatially uniform state* In the discussion immediately 
after Eq. (49), the absence of such a state for sufficiently 
la 0 energies was noted* It was noted some time ago by Onsager 
[46] that, for different reasons, no equilibrium states exist 
for high enough energies. 

Onsager ’s prediction was made in the context of the hydro- 
dynamic formulation of the problem, Eqs. (59) through (61). We 
can understand it in the following way. Equation (3) for the 
partition function can be rewritten as (again omitting factors 
which are not of interest here) 


Z 


■/ 


dX exp (-e/e) 


d£ 0(e) exp (-e/e) 


(115) 


where € is now just the numerical value of the Hamiltonian (59) ; 

2N 

d£ = II (dq^ dpj) is the phase space volume element for all 

the vortices (charges) of both signs, and ft(£) is the structure 
function , or phase space volume per unit £> ft(£) = d$ (£) / d&, 

where $(£) is the total phase space volume with values of H 
less than 6: $(£) = jl*> d-£ f 

For most thermodynamic systems, ft(fi) increases rapidly 
with increasing £ for very large N, and exp(-£/0) decreases 
rapidly for very large N. The competition between the two 
means that the integrand in (115) will be rather sharply peaked 
about some value £ = Indeed, this is nothing but the state 

ment that the fluctuations in energy in the ensemble are small, 
and that most of the systems are sharply concentrated in energy 
near the mean value, which is the only reason statistical mechan- 
ics works at all. Thus we can write 

z = exp ft(e) - e/e) 


(116) 
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where <K£) = ln£2(C). The integrand in (116) can be well ap- 
proximated near £ 0 by exp[tjj(e 0 ) - So./6 + ty"(e. a )(Z ~ e 0 ) 2 /2], 
and is negligibly small away from £o • The integral can be 
performed to give 

-k 

z ss exp t+(e 0 ) - e o /0] 1 t" (e Q ) | /HT (ii7) 

but the theory of the canonical ensemble also shows 


Z 


exp 


*S(£ 0 ) 

K 


fi^/0 + const. 


1 

. 


(118) 


where is the thermodynamic internal energy of the system 
being represented and S(£ 0 ) is the corresponding entropy. The 
maximum occurs at 1/0 - i/j"(£ q ) ■ ^"(C Q )/^(e ). Noting that 
and S are both 0(N) quantities, comparing (118) and (117) gives 
for the entropy 


s(e Q ) = k Urn, n(e 0 ) + const. 


(119) 


plus terms negligible compared to 0(N) , and for the temperature. 


_1 

e 


n'(£ c ) 

*n(e7 


( 120 ) 


Equations (119) and (120) are standard formulas but they have 
some peculiar implications for the system described by Eqs. 
(59)-(61) • For notice that each one of the 2N particles (vor- 
tices) has available to it a maximum phase space volume L 2 |k[. 

All 2N of them have a phase- space volume L 4 n Jk| 2n . This means 
that $(£') , which is a non-decreasing function of £ by definition, 
increases monotonically from 0 to L 4 N|k| 2 N as £ goes f rD m 
to +°°. This further implies- a maximum (at £ = say) in 
ft(£) , where ft'(g) goes from positive to negative* Equation (120) 
then implies that if £ o > £ m > the temperature 0 is negative. 

Landau and Lifshitz [47] show quite generally that a sta- 
tionary thermal equilibrium state is not to be expected for 
negative temperatures, provided the component parts of the sys- 
tem are free to move with respect to each other. That condition 
is fulfilled for the present system, so that no stationary state 
is expected when the value of H exceeds the quantity t , One 
can get a qualitative picture of what the states of tlie system 
are for H - &> £ m by considering the limit of very large 
values of H * These clearly correspond to states in which the 
charges (vortices) of unlike sign as they can get. The 
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macroscopic flow pattern is that of a pair of large counter- 
rotating vortices, spatially separated in different parts of 
the volume. 

Recently, such a flow pattern emerged [37] rather unex- 
pectedly in a numerical solution to the two-dimensional equa- 
tions of incompressible fluid flow at high Reynolds number 
(essentially Eqs. (99) and (100) with a small viscous damping). 
An !l ergodic boundary" separating those initial conditions 
where the large pair of vortices is formed from those where 
they are not was proposed by Deem and Zabusky . A possible ex- 
planation in terms of the "negative temperature" instability 
was offered by Montgomery [48], who also observed that the 
prediction was not totally applicable to the Deem-Zabusky 
simulation in view of the lack of a distinction between "self" 
and "interaction" energies in the continuum representation 
(Onsager’s prediction applies only to a line vortex model, 
and H = G is the interaction energy). Cook and Taylor [49] 
had previously offered an explanation of the Deem— Zabusky 
results in terms of the impossibility of relaxation to the 
thermal equilibrium energy spectrum [essentially Eq, (44)], 
as compatible with conservation of ens trophy > for certain sets 
of initial conditions. [Ens trophy is the volume integral of 
the square of the vorticity, or for the plasma interpretation, 
the square of the charge density. It is one among many con- 
stants of the motion for the ideal incompressible fluid, where 
the integral of any power of the vorticity is also such a 
constant. Due to the non-s quare- in tegr ability of the delta 
function, it is ill-defined for the discrete- charge or dis- 
crete-vortex system] • 

Joyce and Montgomery [50] performed a simulation of the 
discrete vortex situation in which the equations of motion of 
individual vortices were advanced numerically in time. At 
high initial interaction energies (energies » the self 
energies, in the particle-in-cell representation), the Deem- 
Zabusky phenomenon reappeared: formation of a pair of large 

vortices composed of many of the smaller . discrete ones. 

A calculation due to Taylor [51] attempts to determine 
the threshold for the negative temperature instability by 

evaluation of the stationary-phase approximation to lnft(e 0 

The rather remarkable result emerges that, when the infinite 
self-energy is taken into account, £ m = 0. [there is also a 
somewhat confusing attempt to connect the result with the pre- 
vious Cook-Taylor [49] calculation for continuous fluids. We 
are unable to appreciate this argument, and shall remark no 
more on it here] . An improved calculation of £ is due to 
C.E. Seyler , Jr., in Physical Review Letters , 37 515 (1974). 
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The result £ m * 0 would be a remarkable conclusion for the 
following reason. Since the plasma parameter does not enter 
the calculation, we may assume that the result holds down to 
t = 0. But £ - 0 is the "Vlasov" limit, or for the fluid case, 
the limit of ideal, incompressible, inviscid hydrodynamics. 

The spatially-uniform, field-free, quiescent state would be 
by this criterion unstable. The interaction energy per par- 
ticle in the Vlasov limit is zero. Therefore, any perturbation 
on the spatially-uniform state would have the effect of adding 
a positive interaction energy per particle. (The energy 
f dx <E 2 >/8 tt is a positive-definite functional of any charge 
density perturbation) . Thus any perturbation from uniformity 
would put the system on the negative- temperature side of the 
boundary. 

Finally, it may be speculated that the energy at which p 
becomes < 0 according to (49) might lie above £ m . This even- 
tuality, which seems unlikely, would deny the existence in 
principle of a stable spatially-uniform thermodynamic equi- 
librium state of the guiding- center plasma in two dimensions 
for any energy. 

From these various calculations, it is to be concluded 
that the guiding center plasma in two dimensions is an extra- 
ordinarily rich statistical-mechanical system, one about which 
much remains to be discovered. It is also one whose patho- 
logies and singular features should make the investigator wary 
until we have developed a better intuitive understanding of 
the system, than now exists. 


III. GUIDING-CENTER PLASMAS IN THREE DIMENSIONS 

(A) Dynamical Description Canonical Ensemble * and Time- 

Independent Fluctuations 

The two-dimensional guiding- center plasma model developed 
in Sections I and XI is of intrinsic interest for statistical 
physics, but any practical significance it may have depends in 
large part upon finding a way to generalize it to three dimen- 
sions. This program is still in its infancy, and though some 
of the results which have been achieved will now be summarized, 
they will very probably have been considerably enlarged upon 
by the time this material appears in print. 

The natural generalization to three <J_imejjsIons involves 
allowing point charges to move with the cE * B/B 2 (Jrift per- 
pendicularly to a constant uniform magnetic field B, but to 
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move according to Newton's laws in^the parallel direction. If 
the velocity of the jth charge is v j , we decompose it (and 
other vectors) into parallel and perpendicular parts* 



* 

b + v< 




( 121 ) 


wher^ S = B/B is in the z-direction, say, and b * v- =0. 

For Vj^> we have J 

v = cf(x.,t) X i/B 2 (122) 

Ji J 

where E(xj,t) is £he electric field evaluated at tlje instan- 
taneous position x. of the jth charged particle (dx-/dt = v*). 
For the parallel motion, we have 


"J -S 1 - e J 5 ' i( V> 


(123) 


Once again, we use the electrostatic approximation (more 
for convenience than anything- else) and determine E by Poisson* s 
equation, or equation, or equivalently, by its solution 


2(x,t) 



(2 - i .) 

J 

t *+ -4 \ 3 

I* - 


(124) 


where the sum runs over all charges of both signs. (Again 
we take N positive ions and N electrons inside a large volume 
V. Now , V = L 3 ♦ ) 

Equations (121)- (124) form a closed system. The dynamics 
they express can be given a Hamiltonian form; this will guarantee, 
among other things, the existence of a Liouville equation and 
a BBGKY hierarchy, with their attendant possibilities for doing 
a systeroatic^perturbation theory. Letting the jth particle have 
coordinates x* = (x.,y«,z.), the Hamiltonian function can be 
chosen to be 




MAGNETIZED PLASMA MODELS 


483 



where the summations run over all charges of both signs, and 
the canonically conjugate Hamiltonian variables for the per- 
pendicular motion are 



For the parallel motion, the canonical variables are just Zj 
and p z .j ♦ It is straightforward to show that Hamiltonian’s 
equations q^ - 3H/3p-j_, = -9H/3q^ (reproduce (122), while 

zj = 3H/3p z j , p z j = -3H/3zj reproduce (123) and the definition 
°f Vj; = vj z . 

Expressed in terms of the positions and velocities , H 
becomes 


H ■ Z ♦'£ *j v V 2 

i<j J 

where is the two-body Coulomb potential e^ 
The canonical distribution of Gibbs will be 


(127) 

ej/|x ± - xj| . 


,-H/e 


eq.. 


JdX 


J-H/e 


(128) 
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where, if H is expressed in the form (127) , j dX means an in- 
tegral over the 2N • . 4 « 8N coordinates (xj ,yj ,Zj ,Vj z ) , where 
j runs over all charges of both signs* 


That the equilibrium phase space probability density has 
the form that it does [Eq. (128)] means that, in three dimen- 
sions as in two, the thermal equilibrium configuration- space 
probability density (using superscripts to label particle 
species, now) has the form: 






exp[-$(N,N)/9] 

J dx^... dx^ dx^... dxj^ exp 


s exp[-§(N,H)/0]/Q 2K 


(129) 


Equation (129) is the same as for the ordinary plasma in three 
dimensions, with or without an external magnetic field. Here, 
$(N,N) = £ e. e./x. . is the total potential energy of N ions 

i<j J J 

and N electrons. All the reduced probability distributions 



s,+s 
r i e 



dX S p +! ’ ' 


dx N 


expf -f(M,N)/9] 
^2H 


will be the same, too, as will all single-time thermal equi- 
librium ensemble averages. This simplifies matters consider- 
ably, since these quantities are familiar ones and have been 
for several years (see, for a review, reference 13) . They 
are, as before, independent of B. Without proof we shall 
list a number of readily-proved results obtained from 

n 0 , 2 j and ni , i . 

Interest frequently focuses on the case n » 1, where 
now n Q « N/V is the average density (particles/cm 3 ) for both 
species. The Debye length X^ in three dimensions is defined 
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by = 6/8TTn 0 e 2 . For € = (4-rrnaAp) 1 « 1, we have ni )0 
n 0jl = 1 exactly 3 and to 0(e), 


- 11/-41 -41 V 

n 2,0 = 1 + P (*1 “ x 2 > 


ee/-*e -*e>. 

n 0,2 “ 1 + P (X 1 “ X 2 


(130) 


, -+e< 

n i,l = 1 + P (x l " X 2 


ii oe ei ie 

Here, the pair correlation p = p = -p = -p is given 
by 


ee,-K ~ e 
P (*) » - -5“ 


* lAp 


t *4 

I x 


,r* ili'X ee,^, 
= dk e p (k) 


(131a) 


where 


P (k) 


16 t\ j n Q 1 + k 2 


(131b) 


From such expressions, other quantities such as <E(x) E(x r )> 
may be similarly calculated* Thijs , the autocorrelation of 
the electrical charge density p(x) is 


<p(x) p(x')> = 4 e 2 n 2 p ee (x - x') 


+ 2 n o e 2 6(x - x *) 
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„ /-K i6*X - ^ 
S (k) e dk 
P 


where 


Sp(if) - 


2 ,2 ,2 

n o e k 


kir* 1 + k 2 


(132) 


(133) 


Also 


<3(2) 2(2')> 


l E (i?) e lk '< x “ x ) dS 


where 


? ( 2 ) 


from which 

ifi 

8 " 


Un o e ii 


1 + k 


'D 

TTT5 


N) 


_9_ 

2 


dk 


(2n) 3 1 + k 2 ?^ 


(134) 


(135) 


*r 

Equation (135) contains a divergence at large k which, as in 
two dimensions, may be associated with the infinite self- 
energies of the 2N particles 


<2 2 > 


self 1 


8rr 


n e 
o 

2 tt 2 


.2 r 


d & 


2 


at 


"(2 "T k 2 X 2 


(136) 


The 'interaction energy can then be measured by 


<a 2 > - <f 2 > 


self 


8 n 


2 


d£ 


(2n) 3 -k 2 ^(1 + k 2 

( 137 ) 
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so that the ratio of interaction energy to thermal energy is 
approximately 


[<3 g > - <a\ elf ]/8n 

Jn 0 8 


1 

24nn o ^ 


- *§• 038) 


The smallness of this quantity is the basis of the weak- inter- 
action approximation in which the kinetic theory of a plasma 
is usually treated. The limit e 0 for the non- thermal equi- 
librium problem (where n c and 0/m are defined in terms of an 
average number density and average thermal velocity) is usually 
called the "Vlasov limit". The non-equilibrium properties of 
this limit have occupied the attentions of more plasma theo- 
riests than any other single topic. 

To render <E 2 > finite, we may cut off the radial k- 
integration in (135) in the usual way at the inverse distance 
of closest approach of a thermal particle, l% aXi = 0/e 2 . This 
gives <E 2 /8TT>/3n o 0 = 2 tt/ 3, which, it should be noted, is not 
small [unlike (138) J for any 'Value of 0 . The failure to dis- 
tinguish between interaction energy and the spurious selfr- 
energy of point Coulomb charges can lead to serious concep- 
tual errors. 


(B) Transverse Diffusion in Three Dimensions [52] 


The coefficient of transverse spatial diffusion, for the 
three-dimensional guiding- center plasma, c^n be written in the 
sam£ way as (67) , in terms of the part of E that is perpendicular 
to B: 


2 r» 

D = ~ <2 (0) • 2 (T)> dT 

1 B 2 1 i 

J O 


(139) 


Again, the problem^! s on^_ of evaluating the electric field 
auto-correlation <E (0) E (t)>, where E(t) is the electric 
field seen at time^x by a "test" ion which is located at the 
z-axis^a_£ T - 0* E(t) is related to the Eulerian electric 
field E (x , t) by 
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£(x,t) = V ^ 7^| exp i £ • [x -Ic^t)] 

* p 1 x k 

S 3 

* y 

where x(t) is the orbit of the test particle. xj(t) is the 
location of the jtlj charge in the plasma at time t. We have 
toFourier- analyze E(x,t) in a large volume V, rather than an 
infinite volume, for reasons which will become clear as we go 
on* Equation (141) will also be written as 


= £ i, <t) J** , 


(142) 


t 


so that (139) becomes 
■ 00 

n* I 
D 


>. 4 f v jo) ■ * >> = iS,j(T, > ^ 
1 » Jo ^ 


(143 


since <E^ > - 0 for k 1 p — for a spatially uniform ensemble, 

k k 

l 2 


If we once again ignore the correlation between the spatial 
locations of the test ions and the plasma particles [it will 
be small if it is of the same order as (131)], we have a re- 
lation analogous to (77), and 


D = % f V <2 _>(0) • t _ > ('0)<e lk ’ X ( T )> dT. (144) 

x B Jo ^ 2 xk Ak 

On^e again, we can evaluate Dx if we can evaluate 
<exp i k • x(t)>. The same methods will be applicable in 
evaluating 
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16 tt e . e . k . . 

- — p - A -^ -4 <exp ik • [x (T) - x (0)]> 
k 2 V 2 k 2 


(145) 


The program to be^followed closely parallels that for two 
dimensions. For x(t) , we have 


x(t) = 3?(0) + C 


r*. 


Jo 


X t dT' 


(146) 


r T r T# e. b . f(r" ) b 


+ h V|| T + 


dT' dT 


// i 


Jo 


m. 

l 


In (146) , the ’’test 11 particle has been assumed to be an ion 
with charge- to-mass ratio e-L/mj. Equation (146) follows from 
integrating (122) and (123) . x(G) is the ^nitial position and 
the initial velocity of the test ion. x(0) can be set 
equal to zero with no loss of generality, v is a statistically* 
distributed quantity which may be assumed to obey the Maxwell- 
Bo ltzmann distribution. 

Evaluating <exp i k • x(x) is no simple matter, and Eq. 
(146) is considerably more complicated than (78). About the 
best that has Ijeen done so far in evaluating it is the follow- 
ing. Because x(t) is the position of a random variable which 
is initially localized^near x(0) - 0, and because its proba- 
bility distribution P[x(t),t] is expected to spread out with 
time. 


{exp ik * x(t)) 



e itx<T) 


P[x(t),t] 


will damp with increasing T. <exp i k • x(T^_ is the "char- 
acteristic function" of the random variable x(t) , in the lan- 
guage of probability. Roughly speaking, the damping will occur 
for two reasons: motion parallel to the field lines, and 

motion perpendicular to them. Only the latter mechanism was 
operative in the two-dimensional case. The latter is inhibited 



490 


MONTGOMERY 


by increasing B, but the former is not. The damping due to 
the parallel motion will be more extreme than that due to the 
perpendicular motion. Therefore, for very strong magnetic 
fields^, it is useful^to forget the perpendicular drift except 
when k is normal to B. This is tantamount to assuming 



if E? • if f. o 


(147) 


The first of these two expressions can be approximated by the 
assumption which led to Eq . (81), but the second is more com- 
plicated. 

That the statistical properties of the parallel electric 
field seen by a particle can be rather different from the 
properties of the perpendicular field can be seen as follows. 
Even if the parallel electric field obeys some jointly normal 
distribution such as (80) , there is the important difference 
that it can accelerate a particle, whereas the perpendicular 
components cannot. While the perpendicular motion can be visu- 
alized as a limit of small increments in position space, the 
parallel motions are a sum of small increments in velocity * 

space. This means that eventually some of the particles will 
be moving very fast. Such fast-moving particles will, as is 
well known, eventually begin losing energy due to the radia- 
tion of plasma oscillations, and due to the relatively infre- 
quent close collisions with other particles. Therefore, the 
electric field seen in the parallel direction can not in gen- 
eral be well represented by an expression such as (80) in which 
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the probability of the electric field seen by a particle is 
given by a velocity-independent expression. 

The parallel motion is closer in some ways to the 
conventional [31] Langevin approximation of the theory of 
Brownian motion, in one dimension. There, a test particle 
in one dimension obeys the equation 


fjf = -F(v) + A(t) 


(148) 


with a frictional "drag 1 ' coefficient F(v) and a stationary 
random force field A [analogous to our SB # ^(t)]. The 
statistics of A can be given by some such expression as 
(80)* For F(v) « f3v, where 6“ is a constant long com- 
pared to the time necessary for <A(0) A(t)> to go to zero, 
the stochastic solutions to (148) are tractable [31] . 
However, we know that for. the plasma case, F(v) is a 
considerably more complicated nonlinear function of v. 
Methods apparently do not currently exist for the calcu- 
lation of such quantities as <exp[ik v(t ! )dt T ]> when 
(148) governs v(t) , with a more elaborate nonlinear form 
for F(v) . Therefore we must admit an inability to treat 
the second expressions in (147) properly over the whole 
range of T, and seek an approximation which will render 
it tractable, even if the approximation is not entirely 
satisfactory. 

One possible procedure is to neglect the contribu- 
tion of the parallel electric field in the second of 
Eqs * (147). We are treating only the low-frequency, 
long-wavelength parts of the electric field spectrum 
correctly ' anyway. For this reason, one may expect that 
the mechanism of free streaming of particles down the 
field lines will be dominant in destroying the cor- 
relations between the initial electric field seen by 
the particle and that field at a later time. The parallel 
motion of a particle will endure for something of the 
order of a mean free path, and we may imagine that the 
auto-correlation length for the electric field will be < 
this mean free path. If we ignore E , we get for 

t * $ + 0, 
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<exp it • x(T)> k ^ o « <exp(i2 ■ S V|| T )> 


(149) 


dvi. 


Dl. 
1 

2tt0 


exp 


-in. vj 

20 *" 


exp(ikv« T) 


2 2 2 2 2 
a exp(-k ( 0T Hm^) = exp(-k ( V. on T /2), 

where we have assumed the test ions to be distributed with 
a thermal equilibrium velocity distribution* v | on = 6/®i 


We shall make the approximation (149) and explore its 
consequences, both for the test particles and the plasma 
particles. Later we shall remark upon the possibilities 
of relaxing it. For the case of purely perpendicular £ ! s, 
a development parallel to that of 


<exp i 1c 



exp 



(150) 


dT 


Jo 


a, 2 <w • \< t 2 » 


Before substituting (150) and (149) into (144), we 

need expressions for the <E*£(0) • EjJ^(i)>, Making the 

same approximations that led to (149) and (150) on the 
plasma charges in (145) gives 
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<2 JO) • t (r) ) S 

J- k j.k 


I 

ij 


-j £ 2 , 2 

16 rr e , e . k 

g" g "^ “2 < exp ” x.(0)]> . 

k V k 1 3 


f ,2„2 2 ^ 

exp I for k|j ^ 0 


I 

Lid 


i6tt e i 


j^s— exp 


2 ,2 

-c k 


d T 1 dTg • B 4 (t 2 )>] 


. (exp i£ • [x.(0) - x (0)]> for k„ = 0 


(151) 

In (151) , = 0/nu is the square of the thermal speed for 

the species to which charge j belongs- We have also implicitly 
assumed that the jointly normal probability distribution of the 
perpendicular electric field is the same for all the particles. 


Finally, we need <exp i it • [2,(0) - t. (0)]>, which is 
-> 1 J 

J d t 1 dx 2 n 2 (x 1 ,x 2 ) e ±k x i2/v 2 , where n 2 is the probability 
distribution (130) of whatever two particles are involved. 

If. i and j are the same particle, the < > is clearly unity. 
If they are different, i # j, then 
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<exp i k * [x.(0) - 3^(0)]) 




(152) 

sgn e i sgn 


(152) is just (129)-(131), expressed in a finite- volume rep- 
resentation. 

Inserting (152) and (151) into the relation 

(t (0) • t (t)> =V (t (0) . 1 (T)><e lk ' x(T) >(153) 

A 1 ^-'2 ik d 

we have 


<E a (0) * ^(T)> 
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I L 

k n yo 


16 tt e. e . 
sL 

2 2 
k d V 



6^(5). 


2n V 
o 


4 


k 2 + K 


2 flgtle i e j 


2 „2 


-k„ V 


exp 


ion 


I Is 

kii=0 


16 e^. 


v 2 


exp 


/ 2 n 2 r 

' -c k ' 

JL 


V 




■f d \ JI <VV ' V' 


. ,* S 1 S 

°Kr W ~ 2n V , 2 v 2 
o k + K„ 


sgn e. 


(154) 


This is an integral equation for <E^(0) * Ej_(t)>. If we^ha^ 
defined the Fourier transform of the Coulorn^ potential <J> (k) 

lim e 2 /2 7 T 2 (k 2 +X 2 ) more carefully, the 6 (k) terms would not 

. — &r 

X-*0 

have been present in <E_^(0) • E^(t)>, and they will not con- 

k k 

tribute to the summation. Hereafter they are dropped. (This 
is only one of the places where one has to be alert for the 
pathologies that follow from the fact that the Coulomb poten- 
tial is neither square integrable nor absolutely integrable, 
and so has no Fourier transform in a rigorous sense, except 
as such a limit as the one just described.) The can be 
performed in (154) and the result is : 
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<^(0) • ^(T)> 


16 rr 2 e 2 n k 2 ( -k? V? t 2 

S-, r-H «p ■ r 

Vk 2 (k + K 2 ) \ 


H vL -4 4= 

exp s + exp 5 


<T 32n 2 e 2 n o f. 

k ^4) - ■ 


-c 2 k 2 r T r 


/ 4T 1 dT A^l) j. ' 


Equation (155) can be simplified by defining, as in the two- 
dimensional case, 

.§- <2 x (0) . f x (T)> S Q a (T) 


» X < T ) S ~~~2 
X 2 B 


r i dT 2 «l (T 2 ' V 


a 2 r a (t) c 2 q x (t) 


Mi 



magnetized plasma models 

These definitions can be used in (155) to give 


A97 


2Q 1 ( T ) 


1( - 2 2 

16 tt e n 


i-T _v 

(k 2 + 


k 2 /k 2 


k„^0 


+ K") 


exp 


2 „ 2 _2 
ion 


-kfl V - t 


+ exp 


V oe x2 


(157) 


exp 


, 2„2 2 

- k H V ion T 


32 rr 2 e 2 n 


Z 2 1 2 exp [-2k 2 R (t)] 

5 k^ t % 


k ir o 


Defining - 167T e c / B VK^, this becomes, replacing t by t, 


a 2 h (t) 

■L 

at 2 


r- r’ 

/ > , .2 2 + e b / -t f ( k ’ 

1 — 'k 1 + k x S L — ik 


ku=0 


t) 


k„^0 


where f(k,t) is a known function of t: 

2 p p p 

.. k exp -k» V. t /2 

4 — rr% 

k 1 + k . 


(158) 


f ^ V L ^ ] (“ k « V oe ^ 

^exp ^ a J + 6X2 g , 


( 159 ) 
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Equation (158) is the analogue of Eq. (86) for three di- 
mensions* It is considerably more complicated because of the 
last term, however, and so far. has not been solved in the gen- 
eral case. It can be integrated once, noting that dRj/0)/dt = 0, 
to yield: 


dR (t) 
1 

dt 



exp[-2 R i (t)] 

I 72*72 

1 + k , h> 





dT 


(160) 


and we may show how to extract the leading term in powers of 
1/B from Eq* (160) . 

The only place the magnetic field enters (158) or (160) 
is in the denominator of the small quantity E b * We seek the 
leading term, in powers of £ b , of the solution to (158) and 
(160). We shall see presently that for large t, RjU) varies 
proportionally to e* 5 t, which makes the first term on the 
b l 

right of (160) an 0(8^) term. Since the second term on the 

right of (160) is 0(e b ), the result follows that it can be 
neglected to lowest significant order, either in (160) or in 
(158) . 

If we neglect f(k,t) in (158), we get 


a 2 r (t) 

-L 

dt 2 


- Sv 


or equivalently, 
/ \ 2 


i 


L 


k r 0 


dt 


exp[-2 3T R (t)3 

2 2 

1+k A *D 




[1 - exp(-2k R )] 

_ t i (161) 

,2 2v 

V 1 + i V 
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Since Kjl and its first two derivatives are always posi- 
tive for t > 0, last term in (161) becomes negligible as 
t ■+ and 



+ terras of 0( e ^) 


(162) 


(163) 


Equation (163) is the three-dimensional analogue of (89). To 
approximate it by an integral for a large volume, we make the 
replacement [as in (90)]: 



X 

where k min - 2 tt/L = 2tt/V 3 is the lower limit of kjL Integration 
which results from the finite box size. 

Carrying out the integration and introducing the definition 
of (163) becomes 



(164) 


Equation (164) is strikingly similar to (91). It contains 
the factors c0/eB, the square root of the plasma parameter, 
and a volume-dependent factor. Only the volume-dependent fac- 
tor is different: it diverges slowly as L + in two dimensions, 

but it approaches zero slowly in three dimensions. This slow 
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approach to zero puts a constraint on the volume in order that 
the second term of (160) really be negligible compared to the 
first. Replacing the second in (160) by an integral (which 

k 

must be cut off at k = 6/e 2 to avoid an unphysical short- 
max . 

range divergence) the condition that the term involving f(k*t) 
be negligible is that 


B 1 

? » 5 


L/2 it 



4 

(165) . 


This can be regarded as a constraint on either B or L, that 
(164) shall adequately represent the coefficient of transverse 
diffusion. The limits of large B and large V are not inter- 
changeable ♦ At fixed volume, however large, (164) becomes 
accurate as B -* But at fixed B , however large, (164) ceases 
to be accurate as V 00 , because (165) ceases to hold* 

Arguments which led to the inequality (92) on the time 
over which true diffusion- like behavior can be observed like- 
wise complicate the picture for the solutions of (161) . 

(C) Problems Associated with a Convergent Infinite- Volume 

Theory 

The form of Eq. (164) leaves little doubt of its inade- 
quacy if we are interested in passing to the rigorous limit 
of an infinite volume, V ■+ ■». One might at first guess that 
the way to obtain the infinite-volume limit would be to drop 
the first term in (160) and simply write 


D (co) = 2 lim 

1 


a r (t) 
a t 


= 2e b V 


/ 


a k 

(2tt) : 



f(i,r) dT 


( 166 ) 
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However, upon substituting (159) into (166), one finds that 
the integral diverges, and so what superficially appears to 
be^a nice expression with the correct "classical" dependence 
-B is actually useless* 

The reasons for this are not difficult to see, but they 
are not easy to correct. The reasons lie in the separations, 
represented bj Eqs* (147), (149), and (150), of the damping 
of <exp ik • x(t)> into purely perpendicular and purely paral- 
lei contributions.^ The point is that as V increases and the 
allowed values of k become more and more dense, this separa- 
tion will become less and less realistic for vej£y sijall but 
non-zero values of k||* The damping due to the E x B drift 
will compete with the free- streaming damping below a value of 
k |j that can be estimated by setting k 2 ( v| on t 2 = k 2 Dit at a 
time when both are of order unity; i.e., below k|j ^ k \ D x /V ion * 
Similarly the parallel electric field which has been uncere- 
moniously dropped up to now from the second of Eqs. (147) can 
become important in the small kjj region* It has been pointed 
out recently by Vahala [53] that mild assumptions on the sta- 
tistics of the parallel electric field can lead to a term in , 
(147) which has the essential temporal behavior - exp{-k^ t 3 
times a quantity which ^ 0 as e -*■ 0}* This mechanism, essen- 
tially a parallel collision damping, is a j£hird approximately 
equal competitor in the role of damping f(k,r) to zero for 
large T at small kjj. 

It should be apparent that a rather complex limiting pro- 
cess, involving the three small parameters £, £ , and 1/V is 
involved; a simple formula which covers all cases is too much 
to hope for in the near future. But formula (164), supplemented 
by the inequality (165) , may be supposed to be a good measure 
of the regime in which the guiding-center plasma exhibits Bohm- 
like diffusion In thermal equilibrium. What are not yet clear 
are the infinite volume or finite 1/B limits in which "classi- 
cal" diffusion may be expected. The "classical" diffusion for- 
mulas, though very familiar, lack convincing derivations In a 
strongly magnetized plasma. One may expect this part of the 
picture to clear up before long, but jumping to premature con- 
clusions is to be avoided. 

(D) A Systematic Kinetic Theory of the Three-Dimensional 

Guiding- Center Plasma 

Since the system is Hamiltonian [Eqs. (125)-(127) ] , it 
will obey a Liouville equation. From this a BBGKY hierarchy 
can be derived, and provides a possibility of doing a system- 
atic kinetic theory via an expansion in the plasma parameter e. 
This in some ways is similar to the theory for a one-dimensional 
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plasma, but the correspondence is not perfect. We shall adopt 
a compact notation, dispensing with the superscripts to iden- 
tify ‘Charged species , as in the previous two sections* 

For the phase space coordinates of the jth particle, we 
write Xj = (x., y jt z., p zj ) = OyP zj >. where P zj = 

The Liouville equation then reads 



(167) 


where the probability distribution D is a function of all the 
Xj’s for all particles of both species, and is assumed to be 
symmetric under the interchange of like particle coordinates* 
Written out in detail, (167) is 



z-component of the electric field at the location of particle j. 


Summations run over all particles of both signs. 
The reduced probability distributions are 


f 
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J » * 2 S+ 1 


d X 
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(169) 


and integrating (168) yields 



+ v 


5 X. 
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d x 
~s+l 


(170) 
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B 


af, 


s+l 


B X, 
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Here, 


p V’ V*l - j *> 

& 1 "a -*■• 1 

¥i 


9. 


ijS+1 


s+l 1 


/ 1 ' K - vi i 


n o = N / v 


and the £ means to sum over the (s+l) st particle, first 
species 

treating it as an ion and then as an electron. Equation (170) 
is a generalization of (96) , but because of the velocity de- 
pendence, greatly exceeds (96) in dynamical content. 


The M Vlasov approximation 


correlations and writing f «■ 


for the jth species, 

!!i , v J*L . SLi 

£ t z d z m . 


11 amounts to ignoring all the 
s 

II f 1 (i) , as usual. This gives, 
i-1 



Sx| 



( 171 ) 
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where 


3(x,t) = 


-ft I 


“o e 


r 


X - X 


ax' dv' 


species 


(172) 


This is the limit of the Vla&ov-Poisson system in a very strong 
magnetic field. An equivalent limit was treated some time ago 
by Harris [54] and Rosenbluth [55] , although only after the 
solution to the full linearized Vlasov-Poisson system had been 
carried out. 

Standard methods* [13] can be applied to solve the 
linearized version of (171)-(172). That is, we keep only linear 
terms in the departure from the spat tally-uniform, field-free 
state. Linearized, Eq . (171) takes the form 


( 1 ) 

af, 

. ir* 
a t 


(i) 

5f\ 

— 

a z 


e E 


m. A 


8f. 


( 0 ) 




S v 


= 0 


(173) 


where i,e stands for ions and electrons. Poisson^ equation 
(172) takes the form 


2(x,t) 



r 


(± en o ) 


i and e 



» f^^(v ) only are the spatially uniform equilibria we 
x , e i , e z 

are perturbing about, and (x, v , t) are the perturbations. 

l , e z 

Assuming E,f^^ - exp(ik • x) and that ff^ (x, v ,0) ** 

i,e z ^ + 

ik »x ik*x 

g. (v ) e , the solution to (173) is (dropping e ): 
i , e z 


* L. D. Landau, J. Fhye,. (U.S.S.F) 10, 25 (1946). 
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/. \ -Ik v t 

f i,e " g i,e (v z> e * 2 


(175) 




+ li! 5- | dT E (T) e 

m. a v z' ' 

i,e z J 


ik A (u) 


Substituting (175) into (174) gives 


ilt ♦ i?(t) = IkE(t) = 4TTen Q J" d- 


-ik v t | 5f 

; [gi (v z ) - <,(r.» . *■ - * [Vir ♦ T?" 

z z 


(17 6) 


r * -ik v (t-T) 

dr E (T) e z z 


Since E z (t) = (k z /k) E(t), (176) is an integral equation 
of the convolution type for E(t) which can be solved by 
standard techniques [13] « The result is conveniently obtained 
as a Laplace transform 


e(t) 


a+i« 

ds 
2 n i 

cr-i« 


E(s) e S ^ 


( 177 ) 


where s is the complex Laplace transform variable, and 
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E(e) = 


■I 


d v 


z 

is /K 


4 rr e 
Tk 


! tg e (v z ) - g.(v z )] 


(178) 


D(ic,s) 


D(k,s) is the dielectric function, defined by (the usual [13] 
remarks about the contours of s and v z Integration apply) : 




•‘Z. 


2 

E 2 


dv g af|j 0 ) (v.)/av. 


S + 


i k v 
z z 


(179) 


The Zj in (179) is a sum over species (two tjirms) and O) 2 .= 
4un 0 e 2 /m, , ^ = i or e* The vanishing of D(k,s) at ^ = PJ 
- ioj(k) + J y(k) determines the angular frequencies w(k) at 
which the medium can support electrostatic oscillations; the 
corresponding Landau damping decrement is -y(k) . Antecedents 
of Eq. (179) were first derived by Harris [54] and Rosenbluth 
[55]. 


Equation (179) becomes identical with that for an un- 
magnetized plasma upon replacing k z by k, v z by the component 

of v along k, and letting f ^ (v ) be replaced by the velocity 

j ^ 

distribution function of that component of v, It follows th at 
from all the previous detailed work on the solutions of D(k,s) * 
we can immediately infer the character of the possible modes of 
oscillation of the 3-dimensional guiding center plasma. 

In^summary , there are two important branches of the 0) 
versus k curve for Maxwellian distributions , one associated 
with electron oscillations and one associated with ion acoustic 
waves. For k « K^, the electron branch has oj - u)p e k z /k, and 
goes over into ordinary electron plasma oscillations for par- 
allel propagation (k z * k) . The ion acoustic branch is heavily 
damped unless T e » Tj_ , in which case its dispersion relation 


becomes U)/k - (Wpi(k z /k) ) / (k 2 + where Kp e = 47m o e 2 /0 e , 

and this also goes over into it unmagnetized value at 
k z “ k. The damping decrements depend linearly, for the 
case y 2 « oj 2 , on the values of 9fl^/9v z (v z = m/k^) , and are 
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large when these quantities are lajfge, The most important dif- 
ference between the solution of D(k,s) =0 here and in the un- 
magnetized case is perhaps the possibility that undamped electro- 
static oscillations at large angles (k z /k « 1) can exist at 
considerably lower values than the electron plasma frequency: 
ta ~ C0p e k z /k. Since to/k z s (u pe /k, these are still essentially 
undamped as long as L 0 p e /k is greater than an electron thermal 
speed* For sufficiently large wavelength (small k) and large 
angles, these can lie below the ion gyrof requency in a real 
plasma and can render the drift approximation more satisfac- 
tory than it would otherwise be. 

Other kinetic theory problems can be approached through 
(170) . For example, it may be expected that kinetic equations 
for the evolution of non-Maxwell ian one-body distributions, 
such as the Bales cu-Lenard equation, might be extracted from 
an expansion in £ for the spatially uniform state. The strong 
similarity to the one- dimensional case suggests, however, that 
a messy, 0(e: 2 ) , calculation might .be required in order to pre- 
dict relaxation to thermal equilibrium. A calculation of the 
spectral density of the fluctuations in the next sub-section 
points to a possible pathological feature of any such expan- 
sion, and re-emphasizes the very special role played in the 
guiding- center plasma by wave numbers with k^/k « 1. 

(E) Fluctuations for the Three-Dimensional Guiding- Center Plasma 

It is possible, for the spatially-uniform three-dimensional 
guiding-center plasma, to calculate the auto-correlation func- 
tion of the electric field in the Eurlerian representation: 


f -4 -4 

<3{x,t) + r,t + T)> = dit d(i) f e i ^ k,r+<UT ) (180) 


The calculation can be carried out one of two ways: either by 
the Rostoker method of super-posed "dressed 1 ' test particles 
[14] or through the two-time BBGKY hierarchy [13], which can 
be shown to lead to the same results. Both methods are by now 
well known, and we shall use the former since it is simpler. 

The calculation is instructive because it illustrates 
clearly the singular behavior associated with those components 
of k which lie nearly in the perpendicular direction t£ B.^The 
unusual physical behavior associated with k << Kq and k * B = 0 
has already been shown to give the dominant contribution to 
t^e^" anomalous" diffusion coefficient (164) . A calculation of 
S^tkjO)) from the conventional expansion in the strength of the 
interaction sheds additional light on the unusual features of 
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this small, crucial region of k^space. 

"Test particles", in the sense they were introduced into 
kinetic theory by Rostoker [14] and Thompson and Hubbard [56], 
are uncorrelated, freely-streaming Coulomb charges carrying a 
modified potential. The modification of the potential is 
achieved by Fourier transforming it in space an<£ time, then 
dividing the transformed potential 0(k,Q)) by D(k,i(ii). Such 
fluctuation quantities as S^k,Lo) are then obtained by super- 
posing uneowetated particles with this potential. The foun- 
dations of the method are discussed, e.g., in detail in Ref. 13, 
where it is shown to be connected with the more rigorous de- 
rivations proceeding from the two-time BBGKY hierarchy. 

For particles of species j , the exact number density for 
non-interacting charges is: 


n (x,t) 

J 




. /n 

For the 3-dimensional guiding-center plasma, v^ q = b v^, but 
can be anything. The Fourier- transformed charge density 


N 

= l 

2 i=l 


(2lr> 


exp -ik 



+ 



(181) 


will give 



(2 it y 


exp[~iS • 


(3? + v? t) 

10 10 


i£' 


(S» 

v io 


V t 

lo 


i t Sr ~ 1 




( 182 > 
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If the various charges are uncorrelated, the only terms 
which will contribute to the electric field spectrum for large 
t, t f are those with i , so 

<n.0tt) n.(2',t + t)> 

J 1 



(2 n) 3 





(183) 



plus terms which do not contribute to Sg. The velocity inte- 
gral for the 3D guiding- center plasma is a single integration, 
so we can write, for the guiding center plasma 


<nj(x,t) n..(x + r, t + T )> 


(184) 


dm exp[iit 


■4 

r + 


i to ' 


where 




n o,1 

(2? UJ 


(185) 


plus terms which do not contribute to Sg, n Q j = Nj /V = n Q 
for both ions and electrons. Summing (185) over species and 
applying Poisson’s equation, we have the spectral density of 
the electric field fluctuations for ir bare tr (i.e., uncorrela- 
ted and non-interacting) particles: 
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f E (iU 





(186) 


("bare” particles) . 

The canonical recipe now requires division of (186) by 
|D(k,iw) | 2 to give: 


(187) 
i 

I \ I | D(k,iu>) | 2 

The corresponding expression for the three-dimensional ^mag- 
netized plasma is [13,14] 


-Z 


2n .e 2 

QJ J 


k.k 

• 7 


s E (S,».) 


B=0 


„ 2 

r~ • 2 n . e . *♦ ,-» 

■ij " ^ 


_1 

k 


F .(-cu/k) 

... V ■ ■ ■ ■ 

I D(2,iw) l 2 


( 188 ) 


where Fj is the velocity distribution function for the com- 
ponent of velocity along k. 

Equations (187) and (188) are similar, but there is a 
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striking difference. Namely, for U)/k z = const*, k finite, and 
k z 0, (187) diverges; ^hereas (188) does not* In other words, 
£L87^ predicts that < | E (k,o)) [ 2 > will be in no sense small for 
k a B, even though the plasma parameter may be small. 

Superficially more complicated, but similar, divergences 
plague the full three-dimensional plasma in a very strong mag- 
netic field [57]. Ther^ exists at the present time no satis- 
factory calculation of S]r(k,u}) for those components of k nearly 
perpendicular to B. The weak-coupling procedures simply will 
not work. In view of the already displayed importance of these 
components for plasma transport properties, the problem of cal- 
culating them accurately appears to be of the utmost importance. 
Because of the discussion following Eq. (166), it is not reason- 
able to expect this calculation to be easy. 

(F) Summary 

We have outlined a basis for the systematic calculation of 
statistical mechanical quantities for a plasma in which the 
plasma particles E x B drift transverse to a strong dc magnetic 
field, but respond according to Newton's laws in the parallel 
direction. The leading term in the coefficient of transverse 
spatial diffusion for a large but finite plasma has been cal- 
culated, and has been shown to fall off as 0(1/B), even in 
thermal equilibrium. Spectral density calculations for the 
electric field fluctuations based on the more conventional weak- 
coupling expansion are shown to lead to an unphysi^al diver- 
gence for those k vectors nearly perpendicular to B, It is 
exactly these components which are involved in "anomalous" 
transport properties. The conventional expansion in the plasma 
parameter is therefore inadequate in such calculations. 


IV. FINITE GYRO RADIUS EFFECTS AND THERMAL RELAXATION IN TWO 

DIMENSIONS 

(A) Preliminary Considerations; Liouville Equation ; Hierarchy 

Though the guiding- center approximation in two dimensions 
has greatly facilitated calculation of thermal equilibrium 
transport coefficients, certain other processes associated with 
finite values of the plasma parameter apparently lie outside 
its scope. Mot importantly, there is apparently no tendency 
toward thermal equilibrium in the conventional sense. Any un- 
correlated spatially uniform state is an exact solution of (96) : 
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setting f fl s = . n n f ± 0 (i) f 1 (e)gives 3/3t = 0 in (96) 
i' e i=i e=l 9 ’ 

if fi o * fo l s const* = 1* There is no tendency, as there is 
for other Hamiltonian systems that obey a hierarchy, for the 
pair correlations to develop and go to their thermal equilibrium 
values [such as (25)]. Moreover, any state with vanishing trip- 
let correlations and with two-particJ.e distributions which de- 
pend only upon the magnitudes of the separations of the particles 
is also a time independent solution of the equations from (96) 
for which ( Si ,s e ) = (1,1), (2,0), and (0,2). Therefore from 
the dynamical point of view, there is no particular reason to 
prefer (43) over any other spectral density. For example, the 
purely random distribution, with all pair correlations = 0 
and <E 2 > E <E 2 > se ]_f has been studied in some detail numer- 
ically by Taylor and McNamara [28] . 

Clearly what is required in order to study thennalization 
processes in a two-dimensional magnetized plasma is a relaxa- 
tion of the guiding- center approximation. By averaging over 
the "fast" gyromotion, one arrives at a system of dynamical 
equations which, among other things, makes the potential energ; 
a constant of the motion* Since total energy is also conserve 
this makes it impossible for a ratio of potential to total 
energy which is different from that demanded by thermal equi- 
librium to adjust itself. This is an additional proof of the 
impossibility of complete thermal! z at ion within the framework 
of guiding-center theory. 

In this section we shall examine the possibility of pre- 
dicting thermal relaxation for a plasma made up of two-dimens ionc 
charges whose equation of motion is 



instead of (52). To keep the description simple, we shall also 
make the popular approximation of a uniform immobile ion back- 
ground, so that only the negative charges are involved in the 
dynamics : all e^ = — e* The hiouville equation is simple to 

write down, and the BBGKY hierarchy which results from inte- 
grating over the phase space of the last N-s charges is 
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where is given by (15). f s = fgCx^ ... x s v s ;t) is a 
function of the 4s phase space variables necessary to d^scrjbe 
the first s particles, and is normalized so that f f dx. dv. , 

j j t* s t 2S s i i 

dx s dv s = v = L . j 

Equation (189) looks much more like the familiar hier- 
archies that have been investigated than does (96) , but the 
magnetic field term introduces some non-trivlal complications. 


It is realistic to confine our attention to the spatially 
uniform and g^rotr^pi^ ljmit, so that = f x (vj/2,t) only, 
and f 2 = f ^ (x - X 2 » V T* V 2 » t ^ onl y- The s = 1 and s = 2 equa- 
tions become, from (189), 
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Another useful way of writing (191) compactly is 


(& ♦ %) 


+ H 2 j f 2 = L 2 f 3 


(191) 


(192) 


where the linear operator H 2 is the Liouville operator for two 
particles f and L 2 is defined by the right-hand side of (191)* 

(B) Boltzmann's Equation for the Two-Dimensional Strongly- 

Magnetized Plasma 

Equations (190) and (191) are general for any two-body 
potential <f> l2 - Their analogues in three dimensions (usually 
minus the magnetic field) are the starting point for the most 
satisfactory existing derivations of kinetic equations : ap- 

proximate equations of the form 3fj/9t = {a functional of f x 
only} which predict relaxation to a Maxwellian f x * Expansion 
in the density, for example, leads to the Boltzmann equation; 
expansion in the potential leads to the Fokker-Planck equa- 
tion; and expansion in the plasma parameter leads to the 
Balescu-Lenard equation ♦ It is natural to attempt these ex- 
pansions on (190) and (191). 

The perturbation theory is not entirely straightforward 
in the three-dimensional case, and is less so for (190) and 
(191). Normally the complications result from the presence 
of disparate time scales (first articulated by Bogolyubov [2]) 
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in the underlying physical processes* This must be incorpo- 
rated into the mathematics either formally [2,58], or by a 
fortunate insight [59] ♦ The complication of more than one 
time scale is present in (190) -(191) also, but another dif- 
ficulty is more fundamental. The conventional kinetic equa- 
tions have in common that the two-body distribution f evolves 
in terms of and what is essentially a two -particle dynam- 
ical process that goes to completion in a time short compared 
with fi/Bfj/St. For instance, in the derivation of Boltzmann’s 
equation [60,61], the two-body process is a collision in which 
two particles approach each other, collide, and separate in 
well-defined stages during the duration of which f x is sens- 
ibly constant. For two colliding charges in two dimensions, 
the situation is much different in the presence of a strong 
magnetic field, as we shall now show. 

Two equal charges in two dimensions obey the equations 





») 


(193) 


where E l2 - -(2e/£) x 12 /x X2 , 
(193) gives, since E 12 = -S 21 


and X 12 = X 1 ~ X 2 ' 


at 


e 

me 


X t 


Adding Eqs . 


(194) 


where V = (v + v )/2 is the center-of-mass velocity of the 
two charges. Equation (194) is the equation of motion of a 
free particle in a uniform magnetic field and so the solution 
for the center-of-mass motion is immediate. 


Subtracting Eqs. (193) gives 
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a V 2.4 e %£ ft 
d t ~ " m E 12 ” me * ^ 

(195) 

. Uefflg 6 V 4 

s im 2 “ me * 

*12 

where v l2 = v x - v 2 is the relative velocity of the two charges. 
Two constants of the motion of (195) are (as can be readily 
verified) 


+2 

m V 






const. 


mb 





e B 
2c 



const. 


(196) 


Expressing x = (r cos 0, r sin 0) in terms of polar coor- 
dinates and eliminating S, we have 

r 2 + "V(r) ,f = t - const. (197) 

where the "effective potential" is 

■^9 e 2 B 2 r 2 k e 2 

f, V(r)" s — 5— + - — hu? . (198) 

2 m r 8mc X 

The solution to (197) is 

dr 

/ e -" v v(ry ,r 




and since "V(r) ,T has a single minimum and -* 00 asr->Oorr-* 00 
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the r motion is always bounded and periodic, with limits de- 
fined by the intersection of lf V(r) n with £. 

Since 0 = P /mr 2 + eB/2mc, 0 is periodic with the same 
period as r, and w 0 can be written as a constant times t plus 
a periodic function of t with the same period as the r motion 
and with average value zero. The motion may or may not en- 
circle the origin, but there is always a maximum radius r which 
is reached periodically* 

This shows that a significant difference exists between 
the basic two-body Interaction in a magnetized two-dimensional 
plasm^. and a three-dimensional one (or a two-dimensional one 
with B «= 0) * There is no time before which one can confidently 
assume that the spatial correlations of two colliding charges 
vanish. This renders suspect many of the assumptions of the 
BBGKY kinetic theory. Nevertheless, alternative equally- 
plausible assumptions have yet to be suggested, so it is still 
of interest to take the kinetic theory as far as possible in 
parallel to its usual form. 

The easiest expansion to perform on Eqs. (190)- (192) is 
in powers of the density. To lowest order in the density, (192) 
is just the two-particle Liouville equation 

(& + H 2 ) f 2*° . < 199 > 

whose solution is any function of the constants of the motion 
of the two-body problem defined by (193) . If we go ahead and 
ignore the pair correlation at t = 0, f 2 2 f x f x then, the 
solution to (199) 

f 2 ( VVVV T ) = f i[Vi(- T )/2,°] fj_[^g(” T )/2jO] (200) 

where v l (-t) , v 2 (-x) are^th^ s<jlu£ions to (193) that lead to 
two particles, being at x 1 ,v 1 ,x 2 ,v 2 , at time t. Substitution 
of (200) into (190) would lead to secularity for reasons which 
are well-known [11,61,62]. This is avoided by the method of 
Dupree [59] or its more systematic formulation in the method 
of multiple time scales [11,58,61]. 

One notes that an equally accurate solution to (199) is 
obtained by allowing f 1 to have a slow dependence on the time. 

Equation (190) predicts - 0(n Q ), and so a slowly-varying 
f results for low densities. If we write a formal dimensionless 
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expansion parameter £ in front of n Q to remind us that it is 
treated as very small, (200) can to equal accuracy be replaced 
by 

f 2<Vl'VV T ) “ f 1 [\(-T)/2,e t] f 1 tfJ(- T )/2,eT] (201) 

since when substituted into (199) , the only additional terms 
generated are of 0(c) . At the end of the calculation, we can 
then let £ 1. 

The canonical recipe would be to substitute (201) into 
(190) after, passing to the limit T 00 , £T finite. That is, 

one would write 

a (it) 



In the B - 0 case, (202) is nothing but Boltzmanns equation 
[2], and only a few additional manipulations are required [6] 
to bring it into the form given by Boltzmann. For the case 
under discussion, however, it is not so simple. The limit of 
the integrand does not exist; it is an oscillatory function 
of t. The period of oscillation does vary from point to point 
in the two-particle phase space so that the limit of the inte- 
gral may exist even though the limit of the integrand may not. 
But this has not been proved. An additional ad hoc step which 
makes the limit of (202) well defined is that of time averag- 
ing, reasonable in view of the success of the Bogolyubov "method 
of averaging” [63,64] in other situations. It consists of re- 
placing (202) by its time average 
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5 f l (v l/ 2 ' 

a(et) 




3 *12 




(203) 


If the limit in (202) exists, it Is identical to (203). 

We propose (203) [or (202)] as the analogue of Boltzmann’s 
equation for the two-dimensional magnetized plasma [12] ♦ We 
suspect that it obeys an H-theorem and conservation laws, but 
have not been able to prove any of these. 

(C) The Weak- Coup ling Approximation; The Fokker^Planck Equation 

A second expansion of (191) which lies in the mainstream 
of kinetic theory is an expansion in powers of the interaction. 
That is, we shall treat the terms containing a two-body poten- 
tial as of one order higher in a formal expansion parameter £ 
than the terms which do not. Corresponding to this (and glean- 
ing an insight from the equilibrium theory) , we expect the 
correlation functions to get successively smaller in powers 
3f I ~2 

of e as well. Since — ^ = 0(e ), 
f x (l) = f^/2, e 2 t) 


fg(l,2) = t^/2, e 2 t) f x (v 2 /2, e 2 t) 
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+ e P(x 1 ,v 1 ,x 2 ,v 2 ;t, e t) 


f ? (l,2,3) = t^l/2, e 2 t) ^(v 2 ^, e 2 t) f^/2, e 2 t) 


(204) 


+ c[f 1 (v 2 /2, e 2 t) p(x^,v^,XyY^' } t, e 2 t) + cyclic 

permutations] 

+ 0(f) 


and so on. 

Dropping terms of 0(£ 2 ) from (191) after substituting in 
(204) gives 



I.-* 
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f^v 2 ^, ft) ^(Vg/2, e^t) 

(205) 


for the pair correlation P. 

The left-hand side of (205) can be recognized as the total 
time derivative dP/dt computed along the noninteracting trajec- 
tories in the tuo-particle phase space . Integrating (205) from 
0 to t and assuming that no pair correlations exist initially, 
we get 
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p 




3 12 (T) • 


4^ 


^(v 2 / 2, e 2 t) f x (v|/2, e 2 t) 


+( 1 * 2 ) 


(206) 


In (206), the prime means differentiation with respect to v*/2, 
and (1^2) means the same term with coordinates of particles 
1 and 2 interchanged. 

Since x l2 (x), v x (x) , v 2 (t) are siijpl^ periodic functions 
of T, readily expressible in terms of x 1 ,v 1 ,x 2 , and v 2 , the 
integral 



X i2 ( T ) ■ V X (T) 



(207) 


can be done. 1(1) has the physical interpretation of being wT 
times the kinetic energy transferred, to lowest order in per- 
turbation theory, by particle 2 to particle 1. As long as 1(1) 
remains « v*/2, the perturbation theory remains valid. 

The integral (207) Is messy, but can be carried out in 
terms of elementary functions. It is performed in detail in 
Vahala's thesis [12J . Taken over an integral number of gyro- 
periods T q = 27Ttnc/eB, it is 


0 , 


2 2 "l2 ‘ ^ X ^12 ) 

x i2 + a 


> 0 


1(1 )«< 


t 4 tt e 2 V 1 * ( v i2 x ^ 2 2 X 12 ^ X V 12^ 


T jjm 


* X 12 + ' 
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< o 




12 


(208) 
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Here, ft = eB/mc. The content of (208) is roughly that only 
particles initially close enough together contribute to P, and 
this contribution grows secularly with time. 

We must confess lack of a systematic understanding of 
what the proper use of (208) in (190) is. The formal substi- 
tution of (208) into (190) leaves a secularly “growing term 
which is unphysical. It was previously dealt with [12] by the 
somewhat ad hoc procedure of cutting off the integration at ar. 
upper limit corresponding to ml(l) ~ 0. Here we shall use a 
time averaging procedure motivated by "stroboscopic" perturba- 
tion theory [65] and obtain the same answer, up to a numerical 
factor. 


The motivating physical Inequality is that in one gyro- 
period, f 1 should have changed by very little. Therefore, 
(190) will be replaced by 
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* 8*2 


P *12 



(209) 


where the bar over f 2 and P indicates an average over a gyro- 

period, Since in the region where it is non- zero, ttit = 
(27re 2 /£m) • ^ • C^ 12 x b)/v 2 2> and since 


P= -1(1) f'(v*/2, e 2 t ) f x (vp2,% S t) 


( 210 ) 


+ (1 «♦ 2) 


we may substitue (208) and (210) into (209) to get a kinetic 
equation 
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a(e t) 


(^/2, e 2 t) 
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^2 dv 2 


B q>. 


12 


a x- a v. 


^ f i( v i/ 2 > e 2t ) f i< V^' « 2t) 


+ (1 


« 2)^J 


( 211 ) 


with the region of integration to run over all the region de- 
fined by the second inequality in (208). 

The integrals in (211) can be performed [12] and the re- 
sult is 


3 f 1 ( v 2 /2, c 2 t) 

&(e 2 t) 


& v n 


dVg ,( V2 ) • 

( 212 ) 


f 1 ( V l/ 2< ® 2<J ) 


where the dyadic Q is given by 

* 4iT 2 n o e 4 (v 12 X b)(v 12 X b) 

Q " " 2 .2 9 


min 


'12 


,2 4 2 9 

4rr n e v, rt 1 - v, rt v. 
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12 12 12 

—2 


12 


(213) 
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Equation (212) is a standard form in which the Fokker- 
Planck equation has been written before. It is easy to prove 
from it conservation of particles, momentum (zero in this case) 
and kinetic energy. It also obeys an H- theorem, and thus pre- 
dicts an approach to a Maxwellian f^ as t 


L<V 2 ' 


G t) 


~2. 
e t 


m 

2 tt 6 


exp(~m v^/20) 


where 0n o is the initial kinetic energy density. The relaxa- 
tion that occurs is predicted to occur on a time scale t_ 
where 


rel. * 


2 

1 1 

f & o 

rel. 4n o ^ 

That t re ^ shall be long compared to other time scales involved 
is a necessary condition for the applicability of the guiding- 
center model to situations other than thermal equilibrium, as 
done, for example, by Taylor and McNamara [28] for the case of 
a random initial distribution (P ® 0) . 

Up to an arbitrary numerical factor associated with the 
cut-off of 1(1), Eq. (212) is the same as that given by Vahala 
and Montgomery [12]. [See, however, reference [67] .] 

(D) The rr Weak y Long-Range ", or Bales cu-Lenardj Limit 

Note that no long-range cutoff of the potential was re- 
quired to render (213) finite, unlike the corresponding situa- 
tion in the unmagnetized case. Once can thus question the 
necessity for going to the "weak long-range" limit usually 
thought necessary to account properly for the long-range part 
of the Coulomb force. Equation (212) may thus be surmised to 
represent the "Balescu-Lenard limit" as well as the "Fokker- 
Planck limit". Nevertheless, the Balescu-Lenard limit was 
investigated in some detail In Vahala’s thesis [12], P was 
re-calculated by adding in terms 
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3 x 


+ (1 - 2 ) 

on the left-hand side of (205) . The calculation is still man-- 
ageable, though very lengthy. The rather surprising result is 
that following the standard procedures, 3f /3t = 0, through 
terms' of first order In the plasma parameter. This not entirely 
expected vanishing of the Balescu-Lenard collision term has as 
yet no satisfactory physical interpretation, and requires ad- 
ditional Investigation. 

(E) The "Test Particle" Problem; Absence of Dynamic Screening 

on Vlasov Time Scales 

One of the most important physical characteristics of the 
unmagnetized Vlasov plasma in both two and three dimensions is 
the separation of "Vlasov" and "collisional" time scales. These 
can be roughly associated with processes which occur at finite 
rates even when the plasma parameter e goes to zero and those 
which do not. In the former category, we can put plasma oscil- 
lations, Landau damping, and relaxation to the "kinetic stage" 
of Bogolyubov, all of which occur on time scales measured by 
the inverse plasma frequency U)” 1 . In the latter class is relaxa- 
tion to thermal equilibrium, which occurs on time scales - of -1 
in two dimensions, and ~ aT 1 e“ 1 /lne~ 1 in three. p 

Many of the novel features of two- dimensions^ strongly 
magnetized plasmas (and the nearly perpendicular k-vectors in 
three) are associated with the fact that these two time scales 
no longer enjoy a clear separation. There are no relaxation 
processes which occur fo£ components of the fluctuation spectra 
with k perpendicular to B that occur on a time scale independent 
of e, 


This can be illustrated by adapting a very standard kind 
of calculation for three-dimensional plasma [56] to the two- 
dimensional magnetized plasma: the linear response of a Vlasov 

plasma to a "test" charge introduced at t = 0. 

The Vlasov equation obtained from (189) by ignoring all 
correlations is 
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It is soluble, as a linearized initial value problem, by stan- 
dard methods. If the initial conditions are^such that a test 
charge e t is introduced at x Q with velocity v Q = (v Q cos 0 > 
v sin 9) at t =0, the electric field which results for t > 0 
is 



where 



( 217 ) 


£he new symbols appearing in (217) are = |eB/mc|, ^ 
k = (k cos a, k sin a) , and the dielectric function D(k,s) 
is 


D(i?,s) 
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is the Besse l function of its argument and f Q = (2 tt 0) 1 . 
exp (~nrv 2 /20) is the two-dimensional Maxwellian; io pe = 47rn 0 e 2 /£m. 


If the charge is stationary, v =0, and kV « , 

|s|« ^ e> (217) simplifies considerably to t e 

•* ie S 

^.( s ) = — o” 

irk d l (219) 



When we Fourier-Laplace invert (219) we get, for the electric 
field produced by the test charge 




2e t^ - «o> 

i\x- 5 o | 2 


u> 

1 + -? 


all t 


( 220 ) 


which is the vacuum value divided by the low-frequency limit 
of the dielectric function 1 + 

If, however, we had considered the limit 0 in (217) 

we would have gotten instead of (220) , 



for ^p e t » 1 and K D |x - x Q | »‘ 1. 

The differences between (220) and (221) are considerable. 
The latter is Debye shielded and has an effective range 
~ K^ 1 = while the former is not, and has an effective 
range of the order of the size of the system. 


We know, of' course, that the field of the test charge 
will eventually be shielded; the thermal equilibrium theory 
tells us that. The point is that it occurs on time scales 
which are not identifiable with Vlasov^ equation, but with 
processes which slow down to stationary as e 0. 
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(F) Comments on the Problem of Spatial Diffusion for Finite 

Lcumor Radii 

Dawson, Okuda, and Carlile have published a preliminary 
[66] account of a calculation of a two-dimensional coefficient 
of spatial diffusion which endeavors to include finite Lartnor 
radius effects. The essentially new result, for which some 
evidence from numerical simulation is presented, is that the 
Taylor-McNamara result [Eqs. (89) and (91) of this article] 
appears divided by the square root of the dielectric function 

at zero frequency, 

ever, does not asp 
Taylor-McNamara development and its generalizations, and a 
number of the steps are little more than order- of -magnitude 
estimates. These fortunately combine to produce exactly the 
right numerical coefficient to agree with Taylor-McNamara in 
the limit B -+■ 00 . A supplementary derivation based on a two- 
particle encounter model of diffusion is used to interpret 
the data for lower values of B and a "classical" [i.e,, 0(1/B 2 ) 
diffusion coefficient is claimed for these lower values. This 
derivation omits contributions from the very-long-wavelength 
components of the fluctuation spectrum, however, and so con- 
tains no volume dependence. Since it is precisely these com- 
ponents that led to the volume divergence in Eq. (91) and in 
the calculations of Ref, 32, one suspects they may play a sim- 
ilar role in the magnetized f inite-Larmor-radius plasma. This 
"classical" expression therefore probably needs, at best, an 
inequality on the plasma volume in order that it be in fact 
the dominant term. 

One can hope for a clarification of these, and other, 
important unanswered questions connected with the finite gyro- 
radius plasma in the next few years. 


|l + T ^ ie < * eve ^°P ment > h° w “ 

re to the deductive character of the 
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APPENDIX I: DERIVATION OF EQ. (81) 

We are interested in evaluating expressions of the type 

I 3 |4l dt, ... P(g 1 ^ 2 ... exp (iV t m • 

J ^ m=1 


(Al) 


where 


, M 

p ( g l ••• ("5] I <i *,) 


(A2) 




and I) is chosen so that 



P(^ \ ■-) 


1 


(A3) 


The a-Lj are positive definite dyadics which satisfy x . • x > 0 

for any non-zero x. J 


The (E x ,E 2 , ♦ . .Em) can be thought of as a column vector 
whose ith element is say. The I? ± j can be thought of as a 
real, symmetric, positive d^figite matrix whose ijth element 
is j , say. Finally the (Xi , X 2 > t ♦ *X^) can be thought of as 
a column vector whose jth element is f , say. We can then 
write the desired integral I as J 


I - 7] 


(H 


r 


exp 


M 

ij 


A i A e i e J + 1 


1 T. l i i i 


) 


(A4) 


The exponentials in (A4) are scalars, and a unitary trans- 
formation of the £j_'s always exists which will make A^j dia- 
gonal in the new representation. Call the unitary trans- 
formation Uj- j . Thus 
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u i;) £ o 


A ii = Z U ik \n^ u \ 


where 6^ , with > 0, all 1- 

The Jacobian of a unitary transformation is unity , so we 
can write (A4) as 

1 ■ ”1 (3 «$)•" (■ £ a j i t * 1 z *3 e 0 (a6> 

where £' - \ U £. . The integrals in (A6) are now elementar 

j k Jk k 






n* -*f/S 


- H n 7= e 

Now n n (TT^/r^) = 1, and 

i 3 


n exp(-^ 2 /4Ap = ex p(“Z A / 4 } ' 


1 - eX * (“Z A i (A-1) iJ V** 


We also note from (A4) that 



MAGNETIZED PLASMA MODELS 


531 


3 2 I 


all 1-0 
k 


* < e i 


. so 

<«i e j> - <*'V 2 

from (A 8) , and 

I = exp/--!- E V e i V *j) 

V U 

If we now revert to the vector notation, 

M 

■t E 

i,J=l 

To^derive Eq. (81) of the text from (A10) , replace X. by 
( ck X B/B 2 ) At, all for i. 1 

Equation (A2) has the property that integrating it ov^r 
all the Ej_ but one leaves a Gaussian distribution in that Ej_. 
However, the converse is not true: that the electric field 

have a Gaussian distribution at any instant is not a sufficient 
condition that the joint distribution have the form (A2) . 

Considerable use has been made of this theorem in sta- 
tistical physics in recent years. See, e.g., the review article 
by A.J.F. Slegert, In Statistical Mechanics at the Turn of the 
Decade , E.G.D. Cohen, ed.; New York, Marcel Dekker, 1971. 


(A10) 
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